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Preliminary mass transfer studies, studies 
bubble shape and area during formation and 
velocity rise after formation were made for nitro- 
gen bubbles water and two hydrocarbon oils. 

Bubble profiles were found fit empirical 
equation the form where and are polar 
coordinates. 

Quantative studies bubble volume and area 
function time are reported. 

limited correlation for the velocity rise 
chain bubbles viscous oil presented. 


Introduction 


(a) Bubble Shape 

any study mass transfer from bubbles drops, interfacial 
area prime importance. During bubble rise this area 

can evaluated quite readily, but during bubble formation 

nozzle the area varies some manner with time. For 

effective study transfer during formation this area must 

known. 

most studies mass transfer during formation, the 
assumption has been made that the volume the bubble grows 
constantly with time and the surface area computed assuming 
the bubble spherical all times. However, any sequence 
pictures taken during bubble formation reveals that the 
bubble not entirely spherical. Especially during the break- 
off period the bubble has elongated neck which 
attached the nozzle. 

For the calculation either surface area volume 
bubble nozzle, some knowledge the shape necessary. 
For stationary bubble pendant drop, the shape had been 
established theoretically Bashforth and who balanced 
pressure forces inside and outside the surface. The final equation 
could not solved exactly; instead, using numerical method 
Bashforth and Adams prepared tables from which the shape 
could drawn. They also derived equation for the volume 
the bubble above given horizontal plane balancing surface 
tension forces against weight and pressure forces. Tables which 
were quite accurate but only for limited range, were prepared. 

The theory given Bashforth and Adams was used mainly 
calculate boundary tension means pendant ‘arly 
workers, such and Ferguson“) found dif- 
ficulties solving similar equations using the measurements 
they made. Not until Andreas, Hauser, and reviewed 
and simplified the method for finding boundary tension was any 
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degree accuracy attainable. Later solidified thei 
findings and prepared extensive tables. 

The work Lane and and 
Freud and gave good summary bubble 
and volume predicted Bashforth and Adams. Several 
these papers, notably Siemes, discussed solutions the equations, 
all being approximate nature. 

must stressed here that the equations Bashforth and 
Adams were for stationary bubble pendant drop. For 
growing bubble dynamic forces would have considered 
for theoretical solution and, yet, work this nature 
has been done. 

less exact method finding drop volume and shape was 
proposed Null and from direct observation 
the shape. 


(b) Bubble Velocity 

Another important factor mass transfer work the mag- 
nitude the time contact, normally calculated from the 
height rise and the bubble velocity. surveying literature 
one finds any amount data the rise single bubbles 
various liquids. However, little information obtainable 
the rise stream bubbles that could occur study 
mass transfer. 

The rise single air bubbles through different fluids 
widely varying phy sical properties was studied Haberman and 
O’brien and Datta, Napier, and 
and many others. These studies covered the range 
Reynolds numbers from 10~* 3000. For lower Reynolds 
numbers, agreement was found with law 
and the theory put forth and 

The theory Hadamard-Rybezynski 50% in- 
crease the velocity gas bubble liquid circulation 
existed inside the bubble. Bond and proposed 
relationship for the critical radius which transition from 
rigid circulating sphere occurred. Garner and Hammerton 
disgreed with this criterion and presented data for 
numerous systems showing the experimental transition radii. 
However, they did not have enough data for correlation. 

Haberman and Morton studied the rise single air bubbles 
various liquids and also found transition from stagnation 
circulation inside the bubbles. The drag curves the bubbles 
fell between the two limiting curves, i.e. the drag curves 
rigid and fluid spheres. was obtained for the 
transition radii. 

Any velocity data reported for stream bubbles had shown 
increase ‘velocity compared that single bubble. 
Haberman and Morton found that for mineral oil with 
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Figure 1—Theoretical bubble shape. 


velocity and 39% were obtained for bubble diameter 
0.34 cm. with the bubbles 7.7 cm. and 3.2 cm. apart respectively, 


Garner and Hammerton also found increase velocity. 
For white oil (viscosity 198 centistokes 20°C.) velocities 
2.5 times that calculated from Stokes’ law were obtained 
frequencies higher than bubbles per second. 

frequency one bubble per second, Coppock and 
found increase 30% velocity air bubbles 
water for diameter 0.2 cm. the bubble diameter in- 
creased, the difference velocity between single and multiple 
bubbles decreased and became the same diameter 0.6 cm. 
Garner and Hammerton stated that Guyer and 
found increases 50% velocity for carbon dioxide 
bubbles water frequency bubbles per second. 

Other observed increases velocity over 
that single bubbles, but systematic study the effect 
bubble frequency, column size, and seal height bubble velocity 
has been made. any mass transfer work single bubble 
would more difficult analyse than stream bubbles and 
with the latter would imperative know the velocity 
rise order calculate time contact. 


Theoretical Discussion 
(a) Bubble Shape 


good summary bubble pendant drop shape derived 
Bashforth and Adams given For stationary 
bubble attached nozzle and equilibrium the liquid the 


theoretical shape 


where and x/sin are the two principle radii curvature 
point (x, 2), the radius curvature point and 
The gas and liquid densities are and 
respectively, the surface tension and the acceleration due 
gravity The and axes are located shown 
Figure Solutions this equation can found the work 
Bashforth and Adams, and lately extended Fordham. 
The volume, above given horizontal plane was also 
derived and the final equation 


Bashforth and Adams also solved this equation but only over 
limited range. 


Null and Johnson proposed different, less exact, method 
for finding volumes drops. They assumed the bubble 
composed two parts, sphere placed right-truncated 
cone passing through the circumference the nozzle. The 
bubble shape did not lend itself too accurately this treatment. 
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Figure 2—Various bubble shapes for proposed 
Equation 


this work new, empirical, and somewhat simplified 

uation for bubble shape proposed. was found that part 
introducing variable power, place the power the 
equation could used for different bubble shapes Figure 
The resulting equation 


where and are polar coordinates. 


For the surface area and volume Equation (3), from 
revolved about the polar axis. 


sin? 


These equations were solved numerical integration using 
650 computer. Table* and filed with the American 
Documentation Institute No. 6343 permit estimation 
surfaces and volumes revolution; these tables the angles 
are reported degrees rather than radians. 


(b) Bubble Velocity 


For single bubble rising infinite medium high 
viscosity with inertia forces small compared viscous forces, 
Stokes’ law will hold. The resistance the motion the 


bubble will 


where the viscosity, the radius, and the velocity. 
mathematical treatment can found Lamb’s 
buoyant force the bubble against the 
resistance motion results 


where the kinematic viscosity. Equation(8), applicable 
gas bubble, known Stokes’ law, and holds for bubble 


and this paper have been deposited Document No. 6343 

with the ADI Auxiliary Publications Project, Photoduplication Service, 
Library Congress, Washington 25, D.C. copy may obtained 
citing the Document No. and remitting $1.25 for photoprints, 
$1.25 for mm. microfilm. Advance payment required. Make 
cheques or money orders payable to: Chief, Photoduplication Service, 
Library of Congress. 
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Reynolds numbers below the bubble behaves rigid 
sphere with slip the surface. 


The motion fluid spheres was studied independently 
and They found variation 
the resistance motion any circulation existed inside the 
sphere. The derivation also given Lamb. This resistance 


formulated 


where the viscosity the continuous phase and the 
viscosity inside the bubble. Since liquid viscosity much 
greater than that gas, Equation (9) becomes 


and the velocity 
1 2 


This 50% increase compared Stokes’ law, Equation (8). 


Scope This Work 


The purpose this work was threefold: obtain data 
the variation bubble surface area and volume during formation, 
find the actual formation time, and obtain velocity data for 
the systems used. 

Using 16mm. movie film, photographs were taken bubbles 
forming nozzle. From the projected bubble tracing measure- 
ments were taken and the area and volume were calculated 
using equation approximating the bubble shape. Also, from 
these films the actual formation time for various bubble diameters 
was observed. 


Ammonium perchlorate 
drying tube 

Flowmeter 

Thermal conductivity cell 

Nitrogen by-pass 

Ice bath 

n-Hexane 

Control valve for frequency 

Column 

Funnel 

Constant temp. bath, 25°C 

Control panel for T.C. cell 
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Constant temp. water for 

TC. cell, 35°C 


from cylinder 


Figure 3—Schematic diagram apparatus. 
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The velocity rise was measured using two different 
methods which will discussed later more detail. 

The systems studied corresponded the mass transfer runs, 
nitrogen bubble containing n-hexane vapor was passed 
through Paraffinol, heavy mineral oil, and also through light 
mineral oil. Some bubble shapes for air water are included. 


Apparatus 


The apparatus used for photographing bubble shapes during 
formation shown Figure The inlet nitrogen stream 
from cylinder (Linde Company) was passed through mag- 
nesium perchlorate drying tube into the flowmeter hich meas- 
ured the volume flow the pressure the gas. The dry 
nitrogen, after passing through the thermal conductivity cell, 
was directed through saturator containing n-hexane 0°C. 
The saturated stream then formed bubbles the continuous 
phase Paraffinol, heavy mineral oil. the desired bubble 
frequency and diameter, movies were taken using mm. 
Fastax camera (Wollensak Optical Company, Rochester, N.Y.). 
Dupont 930A reversal type film was used speed 1000 
frames per second. The lighting consisted single No. 
photoflood placed directly behind the column, with translucent 
plastic sheet diffuse the light into even background. 


was found that only one nozzle, inside diameter 0.059 
cm., was necessary produce 50% increase bubble diameter 
from about 0.3 0.45 cm. 


The apparatus used for determining bubble rise velocity was 
the same that the mass transfer runs since the data were 
recorded with the transfer data. measure the distance be- 
tween successive bubbles, two techniques were employed, movie 
and still film, and cathetometer (Gaertner Scientific Corpora- 
tion, Chicago, U.S.A.) and stroboscopic light. The frequency 
were obtained from strobotac (Model 631- General 
Radio Company, Cambridge, Mass.). 

was necessary have two reference marks available for 
the film technique. brass rod set into circular thick, 
diameter brass plate was made. The rod had marks 
approximately cm. apart. The complete piece fitted into the 
nozzle assembly. 

Physical properties the materials employed are listed 


Table 


TABLE 


Normal Hexane 


Refractive index, 25°C. 1.3712 
Density, gm./cc. 20°C. 0.65813 
Heavy Mineral Oil 

Gravity °API 28.0 
Specific gravity 60/60°F. 0.8871 
Refractive index, 20°C. 1.4812 
Viscosity, SUS 100°F. 356.9 
Viscosity, SUS 210°F. 52.62 
Viscosity index 
Characterization value 12.35 
Average molecular weight 425 
Composition: Weight per cent 

Aromatics 

Naphthenes 

Paraffins 

Light Oil 

Gravity °API 37.1 
Specific gravity 60/60 °F. 0.8393 
Viscosity, SUS 100°F. 70.3 
Viscosity, SUS 210°F. 36.2 
Viscosity index 111 
Mean average boiling point °F. 705 
Average molecular weight 320 


9 
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Procedure 
(a) Area and Volume 


starting the equipment the flow rate was set approx- 
imately the desired flow and then the continuous phase was 
added the column. This precaution was necessary prevent 
any leakage into the nozzle. After the addition the oil the 
flow rate and frequency were adjusted give the bubble diameter 
needed. small valve just before entry the gas into the 
column served good regulator for the frequency. The 
apparatus was then operated until flow rate and frequency 
fluctuations had disappeared. With feet film the camera, 
movies were taken bubbles forming the nozzle. 


Five bubble diameters from 0.3 0.45 cm. were studied. 


(b) Velocity Rise 

Knowing the frequency and distance between successive 
bubbles, the velocity rise can calculated. This distance 
was determined two ways. 

each bubble moved the column was illuminated 
the light from the strobotac. When the strobotac and bubble 
frequencies were the same, the illuminated bubbles seemed 
remain one place. Actually each bubble moved successively 
from one lighted position another the strobotac flashed. 
Using the cathetometer the distance between these lighted bub- 
bles was measured. 

Any fluctuation frequency either the bubbles the 
light made difficult obtain stationary image. For this 
reason the distance between bubbles was measured photograph- 
ically. Both mm. still and mm. movie film were taken. 


Figure 4—Sequence bubble shapes. System: Paraffinal- 
n-Hexane-Nitrogen; Orifice diameter 0.059 


Equation Final 0.370 em. 


Projecting the film onto screen allowed the distance 
calculated using the two reference marks. 

should pointed out that neither the acceleration nor 
the terminal velocities the bubbles have been measured 
this study. 


Discussion Results 
(a) Bubble Shapes During Formation 


From the mm. film bubble, outlines were traced various 
times during formation, projecting the film onto graph 
Only section film having complete bubble growth and 
release was considered. The time for the various shapes was 
calculated from light exposed the edge the film 
during running operation rate 120 per second. The 
sequence pictures Figure shows the bubble shapes for 
the system, Paraffinol-n-hexane-nitrogen for final bubble diam- 
eter 0.375 cm. 

Since small nozzle was used, the bubble shape was quite 
spherical. the bubble started forming had flattened 
appearance oblate spheroid and this was prevalent with 
all the bubble diameters studied. 

Figure shows bubble shapes for nitrogen bubble forming 
water. The shapes here differed from those Figure 
that the neck the bubble larger and more elongated. Surface 
tension and nozzle size caused these differences. will 
shown later the area variation during formation depends alto- 
gether the bubble shape. 

Figures and values using Equation (3) were also 
plotted show the comparison the actual the predicted 


shape. 


Figure 5—Sequence bubble shapes. System: Water 
Nitrogen; Orifice diameter 0.206 cm.; Equation 
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Figure 6—Sample calculation surface area and volume. 


(b) Sample Calculation Area and Volume 


Referring Figure tracing the bubble outline, the 
ratio was determined. Andreas, Hauser, and 
proposed this same ratio shape factor for pendant drops. 
The values for Equation (3) were calculated for each corre- 
sponding ratio. (Figure 

reference measurement the outside diameter the 
nozzle was used except when the smallest nozzle (d; 0.059 
cm.) was the column. For the latter case greater accuracy 
was obtained using mm. cylinder, mm. high, which 
fitted over the nozzle leaving mm. the tip exposed. 

sample calculation for one bubble shape, nitrogen bubble 
water, Figure was follows: 


value (Figure 2.14 
Actual outside diameter nozzle 0.308 cm. 
Outside diameter nozzle measured 
tracing units 
Maximum bubble width, tracing units 
Actual bubble width 47/42 0.308 0.345 cm. 
Dimensionless bubble width, for 
Equation (3) (Figure for 2.14) 1.715 
Conversion, d/d’, factor change 
dimensionless units actual length 
0.345/1.715 


Angle 
Dimensionless 0.587 
tracing 0.587/1.715 16.1 units 
Value from tracing 16.1 47.9 units 
Value 14.5 units 


0.201 cm. 


Using values for and 2.1 and 17° respectively, the 
area and volume were found from Tables and 


Area 11.37 square units 
11.37 (0.201)? 

Volume 3.744 cubic units 


0.459 cm.? 


the bubble began break away, elongated neck was 
formed which rapidly collapsed until rupture. calculate 
areas and volumes for these shapes, the neck was assumed 
made frusta right cones. 

For the system Paraffinol-n-hexane-nitrogen the odd shapes 
the beginning formation were assumed oblate spheroids 
composed one-half sphere and one-half oblate 
spheroid. 


Variation Bubble Area During Formation 


From each length film for the five different diameters, 
bubble shapes were traced and the area and volume calculated. 


The Canadian Journal Chemical Engineering, August, 1960 


value 


value 


Figure 7—Values for bubble shape. 


For the system Paraffinol-n-hexane-nitrogen, Figures and 
show the volume and area variation respectively. The time 
scale was plotted with the formation time and 
the time starting bubble growth. Since the frequency 
varied slightly this method plotting was preferred. 


None the curves were linear approximate variations 
with time were sought this was important for the mass 
transfer study. The volume varied time fractional power 
indicated the equations and curves drawn Figure 
the particular exponents are empirical. Then, assuming the 
bubble spherical during growth, meant that the area should 
vary with time two-thirds the power used for the volume. 
This proved correct Figure 

However, show that the previous conclusion would not 
hold every case, plots volume and area variation for 
nitrogen bubble water were shown Figures and 11. 
The volume growth was quite linear with time which would 
have suggested that area varied two-thirds power time. 
This was not the area varied quite linearly for most 
the formation time. 

was for this reason that bubble growth was investigated 
thoroughly. For the mass transfer work definite difference 
the amount transfer during formation was predicted depend- 
ing the area variation. 


(d) Formation Time 


From all the film taken bubbles forming nozzle, 
the formation time, was recorded. Figure this time 
was plotted against the flow rate. first was assumed that 
the formation time would vary with frequency but check 
frequencies 151 and 420 proved otherwise. The flow rate 
was the most important variable. 

The points for Figure were obtained three separate 
days and there was vertical deviation shown. Even though 
extreme care was taken cleaning the apparatus, especially the 
nozzle, some contamination must have made the results deviate. 
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Figure 8—Volume variation during formation. System: 
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Figure 12—Formation time versus flow rate. 


However, any particular day the formation time decreased 
with increasing flow rate. the lowest flow rate, i.e. smallest 
bubble size, the valve the inlet line for controlling the fre- 
quency, had closed almost completely. This gave smaller 
formation time, deviating from linear plot through the other 


points. 


(e) Bubble Rise Velocity 


(1) Heavy oil Paraffinol For the system 
hexane-nitrogen, typical curves bubble velocities different 
seal heights and bubble diameters were plotted Figure 13. 
the seal height was increased, the velocity increased, which 
suggested that the continuous phase was moving upward the 
bubble path and was helping the motion the bubbles with 
respect the column. higher seal heights, when the circula- 
tion the column became more developed, the bubble velocity 
remained constant. Actual observation small trapped bubbles 
the continuous phase confirmed the circulation the oil. 

The velocity curves Figure were quite similar and 
interesting find what the velocity would zero seal 
height when theoretically the oil phase had effect the 
bubbles. Several studied the rise single air 
bubbles hydrocarbon oils similar that this work, and 
they found that the bubbles did circulate inside. They assumed 
this circulation since the velocity was that predicted the 
Hadamard-Rybezynski theory. This would seem suggest 
that the intercepts igure should those calculated from 
However, was discovered that the intercepts 
could better represented Stokes’ law shown the 
dashed lines igure 13. From this very useful but limited 
correlation was proposed for this study. 


general curve the difference between actual and Stokes’ 
law velocity was plotted against seal height. later modifica- 
tion changed the seal height the height rise, i.e. the true 
distance the bubble rose. This height rise was calculated 
subtracting from the seal height the distance from the top 
the bubble the nozzle tip immediately after release. Figure 
was prepared from the movie film bubble formation and 
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Figure 13—Bubble velocities. System: Paraffinol-n- 
Hexane-Nitrogen. 


gave the maximum height the bubble above the nozzle 
mouth after release. 

Figure shows the velocity difference from Stokes’ law 
versus the height rise. The majority this data was 
frequency approximately 270 per minute correspond with 
later mass transfer work. Several points were recorded low 
seal heights for the other frequencies needed. Some the mass 
transfer work was frequency 192 per minute and the 
velocity data were estimated from cross-plot velocity dif- 
ference versus frequency with height rise the parameter. 

The best curve through the data 270 per minute was 
calculated least squares analysis assuming polynomial dis- 
tribution. solution was easily obtained from prepared 
program for the 650 computer the Physics Department. 

The curve Figure did not pass through the origin since 
acceleration the bubble took place low seal heights. 

was assumed that circulation still existed the bubbles 
and any agreement with Stokes’ law rather than the Hadamard- 
Rybezy nski modification was superficial. The range Reynolds 
numbers for the heav oil was from 0.6 1.8 which was just 
the transition the curve the drag coefficient versus 
Reynolds number. The drag curve began deviate from Stokes’ 
law this point. Also, and found that 
bubbles behaving like fluid spheres with velocity predicted 
Hadamard-Ry bezy nski, had their velocities decreasing towards 
the solid sphere velocity curve the Reynolds number region 
from 0.1 10. Since the author’s work was this region the 
correlation was probably fortuitous. 

(2) Light oil For the system light oil-n-hexane-nitrogen 
velocity data for bubble diameter 0.34 cm. were plotted 
Figure 16. This curve was essentially the same that 


| | 


Height from tip top bubble 
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Figure 14—Height bubble just before release. 
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Figure 16—Bubble velocity. Light oil. 


Figure 13. However, the velocity zero seal height was cal- 
culated from plot drag coefficient versus Reynolds number 
(Perry’s Handbook 

(3) Uno and studied the rise air 
bubbles through several liquids and found very little wall effect 
the ratio column diameter bubble diameter was ten 
greater. Although the Reynolds number range their work 
was greater (above 20) than this study, since with the largest 
bubble diameter used this work the ratio was greater than 15, 
was assumed that wall effect was present with either 
column used, that the round and the square column 


Time Rise 

For later mass transfer work the contact time for the rise 
period was needed. This was calculated dividing the height 
rise the actual velocity the bubble. From the bubble 
diameter, Stokes’ law velocity, and Figure 15, the contact time 
the bubble for certain seal height could estimated. 

From movie film for bubble diameter 0.40 cm. the actual 
observed contact was recorded and plotted with the cal- 
culated rise. Figure shows the greatest deviation 
high contact This deviation became much less the 
lower umes which the important mass transfer work was done. 
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Figure 15—Bubble velocities different 
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Figure 17—Observed versus calculated time rise. 


Conclusions 


Bubble shapes during formation were predicted quite well 
with variable power. This equation may not 
accurate for static bubbles the theoretical one derived 
Bashforth and Adams but had less complex form. 

many instances bubble volume growth during formation 
had been assumed spherical and linear with time. The surface 
area computed from spherical bubble would then vary the 
two-thirds power time. This assumption would not correct 
every case since surface tension, dynamic forces, and nozzle 
size could combine give very different results. Each study 
bubble formation, especially mass transfer work, should 
include the determination area variation since the amount 
mass transferred will affected. 

limited correlation was found for the velocity rise 
chain bubbles heavy mineral oil. The actual velocity 
was greater than Stokes’ law velocity and varied with seal 
height. The difference between the two each seal height was 
constant the range bubble diameters studied, i.e. 0.33 
0.45 cm. was difficult determine any circulation existed 
inside the bubbles but based the work 
similar systems, circulation seemed very likely. 
The with Stokes’ law was therefore fortuitous. 
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N Re 
a r 
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Greek 
Refe 
(2) 
(3) 
(5) 
4 


(6) 
bubble surface area, cm.? 

radius curvature for theoretical bubble shape 
bubble diameter, cm. 


b 
d = 
bubble width during formation (8) 
measurement bubble shape. (9) 

Reynolds number. 
bubble radius, cm. (13) 
formation time, sec. 
bubble velocity, cm./sec. (14) 
bubble volume, cc. 
Greek (16) 
shape factor. 
angle between polar axis and nozzle edge. (17) 
liquid viscosity, centipoises. 
viscosity inside sphere, centipoises. 
kinematic viscosity, cm.?/sec. (18) 
gas density, gm./cc. (19) 
liquid density, (20) 
radius curvature for theoretical shape. 

radius curvature for theoretical shape. 
surface tension, dynes/cm. (22) 
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Kinetics the Vapor-Phase Oxidation 
Napthalene over Vanadium 
Catalyst’ 


Rates oxidation naphthalene were measured 
for the temperature range 300°C. 335°C. using 
thin layer vanadium oxide-potassium sulfate- 
silica catalyst flow reactor. The amounts 
naphthalene converted were determined analysing 
the products for 1,2-naphtho- 
quinone, and phthalic anhydride. Negligible amounts 
products complete combustion were formed. 

The reaction rate data obtained low conversions 
naphthalene were correlated the following ex- 
pression, based upon model for the reaction 
proposed Hinshelwood: 


From the results experiments which the 
contact time was varied, was suggested that the 
reaction course consisted two parallel, consecutive 
reactions which and are 
intermediates the formation phthalic anhydride. 


the past years numerous experiments have been 
carried out the oxidation naphthalene the vapor- 
phase over contact catalyst. large number these experi- 
ments were concerned with the catalytic properties various 
preparations and the determination conditions for producing 
phthalic anhydride. 

The earliest these studies indicated that the main products 
the reaction were naphthoquinone, phthalic anyhdride, maleic 
anhydride, carbon dioxide and water. Since these can arranged 
the parent compound naphthalene, was usually supposed that 
these products were formed through series consecutive 
reactions. Direct evidence support this view was lacking, 
however, since these early experiments conversions were 
and temperature control the reacting mass inadequate. 

Beginning with the experiments reported 
1952, several investigators have contributed towards 
understanding the kinetics this reaction. Although there 
considerable agreement that the reaction first-order with 
respect the concentration oxygen, the order with re- 
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spect the naphthalene concentration zero according 
while and assumed this order zero and 
that inhibition the reaction products was taking place. Mars 
and van the other hand, treated the reaction 
terms two processes, the first being first-order the naphtha- 
lene concentration and the second first-order the oxygen con- 
centration, with varying steady-state concentration oxygen 
the catalyst surface. Since the conditions used these studies 
varied considerably possible that these differences 
opinion concerning the reaction order with respect the naph- 
thalene concentration are more apparent than real. 


Apparatus 


Figure flow diagram the apparatus used for the 
oxidation experiments. mixture containing known amounts 
oxygen, nitrogen, and naphthalene entered the bottom 
the reactor. These gases passed through the thin bed catalyst 
and then left the reactor through glass head and passed 
air condenser and cold trap. Filters cotton were placed 
between the air condenser and the trap, and near the exit from 
the trap. drop two paraffin oil low volatility was 
usually placed the last filter. 

The carburetors were made from mm. Pyrex tubing with 
porous glass disk near the bottom. 

diagram the reactor shown Figure This was 
constructed from 20-inch length schedule 40, stain- 
less steel pipe located centrally within 12-inch length 6-inch 
steel pipe, and the annular space filled with aluminum. The 
upper end the central pipe was machined take 29/42 
standard glass taper used make connection with the glass 
reactor head. ordinary pipe flange the lower end the 
reactor tube served join this with short piece Pyrex 
pipe and its flange. The gasket between the two flanges was 
Teflon. The outside the 6-inch pipe was covered 
layer mica, over which was placed the electrical resistance 
windings Chromel ribbon. These were then covered 
thin layer Zauer-reisen cement, followed 2-inch layer 
magnesia blocks. 

piece 100-mesh, stainless steel screen was pressed into 
the central pipe serve support for the catalyst. 

The temperature was measured two points calibrated 
Chromel-Alumel thermocouples. The thermocouple inside the 
central tube the reactor was protected sheath made 
from mm. Pyrex tubing. This left narrow annulus between 
the sheath and reactor walls aid removing quickly the 
product gases from the system. During oxidation run the 
lower end the sheath was about inch above the 
catalyst bed. The temperature indicated this thermocouple 
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Figure 1—Flow diagram apparatus for catalytic oxida- 
tion naphthalene. 


was recorded the temperature the catalyst bed. The 
second thermocouple was located near the outside surface 
the central pipe means hole drilled through the aluminum 
block. Bristol model 560 Dynamaster was used record 
the temperature indicated each thermocouple and regulate 
the power input the electrical resistance windings the 
outside the reactor. The sensitive element this control 
was the thermocouple the aluminum block. For most the 
oxidation runs the two thermocouples gave the same temperature 
reading; times, however, the temperature the aluminum 
block was about degree higher. 


Operating Procedure 

The operating procedure involved, the first place, pre- 
liminary oxidation run 1-2 hours. The purpose this was 
permit the system, particularly the catalyst, attain steady- 
state condition. run proper was then started replacing 
auxiliary condenser weighed air condenser and switching 
the primary gas stream from the unweighed naphthalene car- 
burettor. The time for these runs varied from three-quarters 
hour two hours. the end day’s operation nitrogen 
was passed through the reactor for least minutes before 
the current the electrical heaters was turned off. 


Materials Employed 


Naphthalene was B.D.H. Certified chemical. Melting point 

Benzene was Fisher Certified Chemical, thiophene-free. 

Oxygen and nitrogen were taken from cylinders containing 
the compressed gases. They were dried over solid potassium 
hydroxide before being passed the reactor. 

1,4-Naphthoquinone (Fisher) was recrystallized from ether 
suggested (melting point 125°C. corrected). 

1,2-Naphthoquinone was prepared from Orange-2 using the 
procedure outlined The material did not give 
sharp melting point but softened and turned black about 
consisted tiny needles colored golden orange. 


The catalyst was obtained from Imperial Smelting Corpora- 
tion Limited, Avonmouth, England. was taken from large 
batch catalyst which the mean analysis was: 

Loss ignition 12.5 

The remainder mainly silica. 
All analyses are expressed percentage the catalyst ignited 
800°C., except the figure which expressed the 
sample dried 110°C. 

All the naphthalene conversion rates reported this paper 
were carried out using single sample this catalyst weighing 
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Figure 2—Reactor for catalytic oxidation naphthalene. 


1.130 grams. The size the particles this sample was kept 
uniform spreading out the crushed catalyst graph paper 
and picking out hand only those pieces which fitted square 
0.05 inches side. 


Identification Reaction Products and 
Methods Analysis 

Phthalic anhydride was qualitatively identified the reaction 
products fusion with resorcinol and formation the dye 
was determined quantitatively titration 
with 0.02 NaOH using phenolphthalein indicator. 
direct titration was found unsatisfactory due the end- 
point being obscured the small solubility naphthoquinone 
the aqueous solution. Therefore, the product solution 
benzene was first extracted with 0.02 NaOH, the aqueous 
layer acidified with 0.02 HCI, this solution extracted once 
with benzene remove naphthoquinone, the aqueous layer 
boiled remove carbon dioxide, and final adjustment the 
end-point made after cooling the solution room 

1,4-Naphthoquinone was identified the reaction products 
measuring the absorption spectrum benzene solution 
the products after this solution had been extracted several times 
with 0.02 NaOH. Beckman quartz spectrophotometer 
was used make the measurements. The position the peaks 
the region 325 600 my, well the relative extinction 
these peaks, agreed closely with the measurements 
reported Nagakura and and determined with 
sample pure 1,4-naphthoquinone benzene. 

During the period when preliminary runs were being carried 
out was noticed that the products contained orange-colored 
substance. The only colored substance expected from informa- 
tion the literature was 1,4-naphthoquinone, which colored 
bright, canary yellow. was also noticed that the products 
stood open the atmosphere, color change dark green 
occurred 2-3 days time. number experiments were 
then made attempt discover the nature this unknown 
substance. one instance, the total naphthoquinone content 
was determined ceric sulfate titration (method outlined 
below) and Klett-Summerson colorimeter calibrated for 
1,4-naphthoquinone. The value determined colorimetrically was 
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several times larger than that found titration. This indicated 
that substance was present which has extinction coefficient 
greater than that 1,4-naphthoquinone the wavelengths used 
the colorimeter (band maximum 420 my). was postulated 
that the unknown substance was 1,2-naphthoquinone. test 
this, the products from runs were analysed for total 
naphthoquinone titration with ceric sulphate solution and 
extinction measurements using extinction determined 
from pure samples 1,4- and 1,2-naphthoquinone. 
these runs the two methods gave the same naphthoquinone 
content within experimental error. another series tests 
synthetic mixtures were prepared containing, one case, naph- 
thalene, phthalic anhydride, -naphthoquinone, and 
second case, 2-naphthoquinone well the other constituents. 

These were let stand open the laboratory. Only 
the mixture containing 1,2-naphthoquinone developed green 
color after 2-3 days. 

The chemical method used analyse for naphthoquinone 
was that previously outlined Rosin ml. benzene 
solution the reaction products was treated with ml. gla- 
cial acetic acid and 0.5 gm. zine dust. When the solution 
was colorless, the mixture was diluted with water, filtered, 
and titrated with standard ceric sulfate solution using o-phenan- 
throline indicator. 

For the extinction measurements Beckman quartz 
spectrophotometer with matched cells silica was employed. 
Solutions 1,4- and 1,2-naphthoquinone, respectively, were 
shown follow Beer’s law wavelengths 340 and 540 mu. 

The product gases from runs the temperature range 
312° 352°C. were analyzed for carbon dioxide determining 
the increase weight absorber containing Carboxite. 
all but two these runs the amounts carbon dioxide were 
found equivalent, within experimental error, the amounts 
anhydride formed. was therefore concluded that products 
corresponding oxidation beyond the phthalic anhydride stage 
were not formed under the conditions temperature and contact 
time employed. 

The moles 1,4- and phthalic 
anhydride the product gave, therefore, the moles naph- 
thalene oxidized during the run. This value, together with the 
number grams cataly the reactor and the length the 
run seconds, was used calculate the naphthalene conversion 
rate moles naphthalene oxidized per gram catalyst per 
second. 
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Figure 3—Rate oxidation naphthalene function 
the average concentration naphthalene the bed 


2.08 moles per liter; 4.64 moles per liter; 
6.86 moles per liter. 
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Experiments with Catalyst the Reactor 


number runs were carried out under conditions similar 
those during the catalytic runs but with catalyst the 
reactor. The products obtained contained acidic substances 
but were colored faintly yellow. Extinction measurements were 
carried out 340 solutions these products benzene 
and the results interpreted terms 1,4-naphthoquinone con- 
tent. Comparison these amounts 1,4-naphthoquinone with 
those formed during the catalytic experiments indicated that the 
non-catalyst formation naphthoquinone one percent less 
that formed during the catalytic runs. correction the 
measured catalytic rates was therefore made. 


Oxidation Runs Low Naphthalene Conversions 


this series runs the results were plotted naphthalene 
conversion rate versus the arithmetic average the naphthalene 
concentration entering and leaving the catalyst bed, with oxygen 
concentration parameter. Figure the results are shown 
335°C. These curves show that the conversion rate dependent 
upon the naphthalene concentration low values the naph- 
thalene concentration and high values the oxygen concentra- 
tion. high concentrations naphthalene, however, the 
conversion rate nearly independent the concentration 
naphthalene. 

Although oxidation the naphthalene most probably occurs 
near the catalyst surface, possible that the rate 
adsorption one reactant the same order magnitude 
the rate the chemical reaction. the following treatment, 
that oxygen only adsorbed the catalyst surface, that the 
rate desorption oxygen negligible, and that reaction 
occurs when naphthalene the gas strikes adsorbed 
oxygen, then can expected that steady-state will 
established wherein the rate removal oxygen chemical 
reaction equal the rate adsorption oxygen. Thus: 


the number moles oxygen consumed for each mole 
naphthalene reacting, then: 


RAM OF CATALYST PER SECOND X 10° 
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Figure 4—Rate oxidation naphthalene function 

the average concentration naphthalene the bed 

catalyst. Temperature 335°C. Oxygen concentration: 

3.00 moles per liter. 6.61 moles per liter; 
8.81 moles per liter. 


From 


For con 


When 


correlat 


Using 
are 
that 
concent 
given 
values 

Add 
concent 
tially 


Since 
calculat 

The 
constan 
imately 
This 
concent 

The 
plot 
tration 
was 
rate 
concent 
conclud 

Sinc 
Lang 
surface 
tried, 


| 
v 
4 v 
6 Oo 
| 
— 
u 
4 
a 
4 4 | 
: 
| 
312 
| 
10. 
‘ 2 


(1/N) (rate reaction oxygen) 


From and (2), the rate oxidation naphthalene written: 


When inverted, the equation easy use for purposes 


correlating the data. 


Using Equation (5) and the method least squares, the most 
probable value for constants 1/k, and respectively, was 
calculated from the experimental data. Using these values 
Equation (4) the vs. curve was calculated. These curves 
are given the full Figures and can seen 
that reasonable fit the data obtained constant oxygen 
concentration. for the constants Equation (4) are 
given Table this table, the average the individual 
values for each series runs constant oxygen concentration. 

Additional values were obtained from runs high 
concentrations naphthalene where the reaction rate essen- 
tially independent the naphthalene concentration. For such 
conditions, Equation (3) predicts: 


Since the values and are all measured, readily 
calculated. Particulars concerning these runs are listed Table 

The specific velocity constants, and k,, are expected 
constant given temperature. the case this approx- 
imately so. However, considerable variation was found for 
This variation does not appear due change the 
concentration oxygen. 

The oxygen dependence the rate was checked 335°C. 
plotting naphthalene conversion rate versus oxygen concen- 
tration for series runs which the naphthalene concentration 
was greater than moles per liter. this range the 
rate oxidation naphthalene nearly independent the 
concentration naphthalene. The curve obtained (see Figure 
essentially linear, extrapolating through the 
concluded, therefore, that the naphthalene conversion rate 
the oxygen concentration. 

Since the shape the curves Figures and similar 
Langmuir-type adsorption isotherms, may argued that 
method correlating the data may possible based 
surface reaction between adsorbed reactants. This was indeed 
tried, and, although slightly better correlation the data was 
obtained, interpretation the constants the resulting equation 


CONSTANTS EQUATIONS (3) AND (5) 


‘TABLE 


681 


2.51 


proved difficult. The expression obtained was more 
complicated than (3). for this reason, together with the 
fact that there considerable scatter the data available 
the present time, that the simpler Hinshelwood approach 
preferred. 


Naphthalene Conversion Longer Contact Times 


additional series runs 335°C. the concentration 
oxygen and naphthalene was kept approximately constant 
while the contact time was varied changing the gas flow 
rate. The results obtained are presented graphically Figure 
the figure the reciprocal the gas flow rate plotted 
the abscissae since this directly proportional the contact 
time for fixed amount catalyst. seen that the concentra- 
tion 1,4-naphthoquinone and phthalic anhydride increase 
almost linearly over the range flow rates used, while, over 
this same range, decrease the concentration 1,2-naph- 
thoquinone occurs. 

interpret these results terms the reaction course, 
supposed, first all, that 1,4- and 1,2-naphthoquinone are 
formed through two separate reactions occurring simultaneously. 
This reasonable since that either quinone can 
formed from the other under the conditions the oxidation 
experiments. the second place, reaction leading the re- 
moval 1,2-naphthoquinone involved since the concentration 
this quinone decreases for increase the contact time 
(see Figure 6). Concerning the formation phthalic anhydride 
three possible sources are considered: 

(1) Phthalic anhydride formed process parallel 

the formation 1,4- and 

(2) Phthalic anhydride formed from 1,4-naphthoquinone. 

(3) Phthalic anhy dride formed from 2-naphthoquinone. 
With regard point (1) that third simultaneous reaction 
gives rise the phthalic anhydride, plot the moles 
anhy dride per mole 1,4 naphthoquinone the product versus 
the ‘reciprocal the gas flow rate should horizontal, unless 
the two reactions are different order with respect the 
naphthalene concentration. Such plot given Figure 
from which seen that this ratio increases with the contact 
time. Therefore, either phthalic anhydride not formed solely 
through third simultaneous reaction, else the two reactions 
are different order with respect the naphthalene concentra- 
tion. the author’s opinion that such difference order 
unlikely. Concerning point (2) that 1,4-naphthoquinone 
intermediate the formation phthalic anhydride, the slope 
the curve Figure should decrease the contact time 
increased. Clearly significant effect this kind indicated. 


phthalic anhydride, suggested point quite possible, 

however. Figure shows that the concentration this quinone 

decreases with contact time expected for intermediate 

two consecutive reactions. therefore concluded that the 

data represented Figures and are consistent with 

scheme involving two parallel, consecutive reactions: 


phthalic 


naphthalene 
anhydride 


TABLE 


REACTION RATES AT HIGH CONCENTRATIONS OF NAPHTHALENE 


t 
_ 
w 
tw tw 


For constant oxygen concentration, Equation (3) becomes: 
(4) 
| | | | | | j 
| | | ; 
312 4.64 242 274 4.14 1.94 2.47 
312 2.32 3.08 2.60 1.96 
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Figure 5—Rate oxidation naphthalene function 
oxygen concentration for concentrations naphthalene 
greater than moles per liter. Temperature 335°C. 
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Figure 7—Relative amounts phthalic anhydride and 

1,4-naphthoquinone formed the catalyst bed after dif- 

ferent contact times. Oxygen concentration: 6.6 
moles per liter; 4.4 moles per liter. 


Comparison with Work Reported Other 
Investigators 


concluded that the rate naphthalene oxida- 
tion zero-order the concentration naphthalene. The 
present study shows this only for high naphthalene concen- 
trations. his plots naphthoquinone formation, Calderbank 
found apparent conversion W/F equal zero. 
concluded, therefore, that part the naphthoquinone formed 
homogeneous reaction. The experiments carried out 
the present author which catalyst was present the 
reactor (discussed above) have shown that the extent this 
‘so-called’ homogeneous reaction much too small account 
for the size Calderbank’s intercept. possible that the 
reaction products obtained Calderbank contained 
1,2-naphthoquinone which was not taken into account his 
method analysis. Since used method 
determine naphthoquinone, likely that his naphthoquinone 
analysis too high due the greater extinction 
and that this the real reason why his 
plots indicated intercept for equal zero. 
also concluded that 1,4- naphthoquinone anhydride 
are formed simultaneous reactions when the conversion 
naphthalene small. The present study agrees with this, with 
the modification that 1,2-naphthoquinone intermediate 
the formation phthalic anhydride from naphthalene. 
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Figure 6—Variation with contact time the concentration 

1,4-naphthoquinone (), phthalic anhydride and 

1,2-naphthoquinone the gases leaving the bed 

catalyst. Temperature 335°C. Oxygen concentration 6.6 
moles per liter. 


Mars and van measured the rate naphthalene 
conversion fluidized bed catalyst. With this technique 
these authors were able realize small naphthalene concentra- 
tions the bed catalyst. They observed that the reaction 
rate was dependent upon the concentration naphthalene, 
was also observed during the present study. These authors 
correlated their data using kinetic equation identical form 
with (3), based upon model the reaction involving 
alternate reduction and oxidation the catalyst surface. 

loffe and used integral reactor their study 
this reaction. They found the reaction rate approx- 
imately first-order the concentration naphthalene temper- 
atures below 350°C. but that the specific rate constant this 
basis decreased value for increase the concentration 
naphthalene the feed gases. This anomalous behavior the 
first-order constant understandable the light the observa- 
tions noted during the present study. Since the naphthalene 
concentrations used loffe and were relatively 
high probable that their reaction rate was more nearly 
zero-order than first-order the naphthalene concentration. 
increase the concentration naphthalene then expected 
result decrease the value the specific rate constant 
calculated the reaction were first-order the naphthalene 
concentration. 

the naphthalene reaction, loffe and Sherman concluded that 
scheme parallel and consecutive reactions was invloved. The 
study agreement with this, with the modification 
that 1,2-naphthoquinone intermediate the formation 
phthalic anhydride from naphthalene. 

and reported the oxidation naph- 
thalene first-order the concentration naphthalene. 
The present study agrees with this since these authors used low 
concentrations naphthalene and high temperatures. 


Conclusions 


Rate data are reported for the oxidation naphthalene 
under conditions resulting low conversions. This has been 
correlated the basis model proposed Hinshelwood 
that steady-state concentration oxygen established 
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the surface. The resulting equation has also been 
shown consistent with the data the literature. 

shown that 1,2-naphthoquinone formed small 
amounts during the oxidation naphthalene. This substance 
has not been reported previous investigators the kinetics 
this reaction. 

scheme parallel and consecutive reactions has been 
suggested for the course the chemical reaction. this 
scheme 1,4- and 1,2-naphthoquinone are intermediates the 
formation phthalic anhydride. 
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Nomenclature 


gas concentration, moles per liter. 

gas flow rate reactor, liters per minute. 

specific rate adsorption oxygen. 

specific rate reaction naphthalene. 

specific rate reaction adsorbed oxygen. 

equal C,). 

moles oxygen consumed per mole naphthalene 
converted. 

reaction rate, moles per gram catalyst per second. 
gas constant. 

fraction the catalyst surface covered with oxygen. 
temperature. 
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Subscripts 


n 


naphthalene. 
oxygen. 
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Oxidation Ethylene 
Using Silver Catalyst Coated Strips 


Fluidized bed tests with silver catalysts used 
the oxidation ethylene ethylene oxide showed 
fluidization instabilities due catalyst agglomeration. 
One method avoiding the agglomeration problem, 
yet utilizing the excellent heat transfer property 
fluidized bed, the immersion fixed catalyst 
coated strips inert fluidized bed. Small scale 
tests with this method, have shown that excellent 
temperature control can achieved without serious 
adverse effects catalyst performance. Kinetic data 
are presented for limited range flow conditions, 
and are correlated empirical rate equation 
the form 


silver catalyzed oxidation ethylene leads the main 
reactions 


Large amounts heat are evolved this oxidation, particularly 
reaction (2). Close catalyst temperature control essential 
order achieve maxiumum yield ethylene oxide, and 
maintain catalyst stability. One major requirement therefore 
any practical reactor design suitable for these reactions 
adequate heat transfer rate. 

The fluidized bed reactor would appear applicable because 
its excellent gas-solid contact and its inherently high transfer 
rates. Its application has been hindered part however 
catalyst agglomeration which cause fluidization insta- 
bility. Packed tube reactors avoid agglomeration effects but 
have the disadv antage relatively low heat transfer rates. 

combination fluidized and fixed catalyst bed may 
however acceptable basis for reactor design. Catalyst 
coated strips can prepared from which are 
mechanically strong, and which can withstand the attrition 
effects fluidized bed. bed material such glass beads 
inert this reaction, and will serve excellent heat transfer 
medium. This combination should have some the advantages 
fluidized bed such high heat transfer rates, and easy 
catalyst replacement. would have similar limitations 
space velocity range, and would perhaps less efficient gas 
iManuscript received February 4; accepted April 2, 1960. 
2Tokyo Institute of Technology, Ookayama, Meguro-Ku, Tokyo, Japan. 
3Division of Applied Chemistry, National Research Council, Ottawa, Ont. 
Contribution from Division of Applied Chemistry, National Research 
Council, Ottawa, Ont. Issued N.R.C. No. 5685. Based paper 


presented to C.1.C. Chemical Engineering Conference, November 9-11, 
1959. 
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catalyst contact. Because the novelty the combination, 
and its potential application the catalytic oxidation ethylene, 
performance data were obtained small scale model reactor 
under variety process conditions. 


Experimental 


The reactor tube shown schematically Figure 
made from in. I.D. in. stainless steel tube; with 
steel ball the bottom the conical bottom act check 
valve. electrically heated, and the fluidized bed within 
kept constant temperature means temperature 
controller. Bed temperatures were recorded from 
couples which passed through the wall and into the bed. The 
temperature difference between top and bottom bed thermo- 
couples was about 2°C. Catalyst strip temperatures were 
not measured; that temperature data this report are based 
the fluidized bed temperature. 

The catalyst sprayed stainless steel strips which are 
attached centrally located supporting structure, and entirely 
immersed the fluidized bed. The strips are arranged four 
stages, each stage containing four strips spaced 90°C. apart, 
and rotated 45° with respect one below it. The dimensions 
strip are 0.63 2.56 .033 in. These strips were rough- 
ened grit blasting, and passivated 50% nitric acid 
solution. The catalyst alloy (8.25% 91.7% Ag) was put 
onto each side thickness 0.02 in. flame spray method. 
Total weight alloy the sixteen strips was 37.6 gm. 

The activated catalyst was obtained steaming 320° for 
hours, followed boiling 20% acetic acid 
washing. The weight activated catalyst remaining was 33.2 
gm. Its residual calcium content was about 0.1% calcium; and 
its initial specific area was about 0.5 sq. m./gm. 
activation the catalyst was conditioned operating the reactor 
for about two Its activity-time pattern during this 
period was similar that reported and 
except that longer period was required 
reach constant level activity. The data reported this 
paper are based performance after the two week conditioning 
period. 

Glass beads, 150 200 mesh, were used the bed 
material. Other materials such artificial graphite and silica 
were also graphite caused minor catalyst 
poisoning, and formed dust. Silica gel dusted readily, and also 
cataly zed the isomerization ethylene oxide acetaldehyde. 
This isomerization could eliminated however treating the 
silica gel with mole sodium carbonate. Glass beads were 
superior both these bed materials, since poisoning occurred 
and elutriation was negligible. 

Erosion the cataly surface the fluidized glass bead 
bed was not significant variable under the flow conditions 
tested. After the initial conditioning period catalyst activity 
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Figure 1—Mounting catalyst strips reactor. 


remained substantially constant over several weeks continuous 
operation, which would not have occurred appreciable amounts 
catalyst were removed erosion. Some erosion may have 
occurred but could not detected weighing the strips since 
glass beads had become embedded the larger holes the 
surface. Visual examination strips the end the experi- 
ments did not reveal obvious wear. 

The air was supplied the reactor from Nash Hytor 
compressor which distilled water was used the sealing 
fluid. Oxygen, nitrogen and ethylene were supplied from 
cylinders without further purification. Ethylene was obtained 
from Ohio Chemical Co. Limited, and was 99.5% purity. All 
feed gases the reactor were metered with rotameters operat- 

Exit gases from the reactor passed through clone, 
pressure reducing valve, then cooler, and finally water 
bubbler, before being vented exhaust system. sample 
exit gas was taken from the line immediately ahead the 
bubbler, and analysed for and CO, gas 
chromatography Accuracy analysis for individual 
components was 2%. experiments where low conversions 
were obtained, Drierite and Ascarite traps were placed the 
inlet air line reduce the possible error product analysis 
caused the presence and the inlet air. 


Correlation Conversion Data 


Conversion data function W/F ratios were obtained 
constant fluidized bed temperature and constant air 
ethylene ratios. The ratio W/F was varied holding constant 
and varying i.e. mounting two, three four sets 
catalyst strips series within the fluidized bed. The flow rate 
was set about three times the minimum required for fluid- 
ization. alternate method changing W/F holding 
constant and The range gas flow limited 
the range between minimum fluidization and the onset severe 
slugging. When slugging occurred, large gas bubbles uncovered 
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the catalyst causing poor heat transfer, and uncontrolled 
catalyst temperature. Within the gas flow range good 
fluidization, adequate temperature control was obtained. Conver- 
sion data the same ratios were not significantly different 
between the two methods changing ratios. 


Ethylene Dependence 


Conversion data based variable and constant are 
Table where the fraction ethylene feed which 
has been oxidized all products, and selectivity, the 
fraction the total ethlene consumed which has been oxidized 
ethylene oxide. Plots the data several bed temperatures 
with air-ethylene ratio 15:1 are given Figure Similar 
plots were made with air-ethylene ratio 20:1 and 35:1. 
Initial rates were obtained drawing tangents best fitting 
curves W/F and these slopes were replotted log 
against log Linear plots were obtained each temperature 
having common slope 0.5, that the initial rate data may 
represented the empirical Equation (3). 


was assumed that the oxygen concentration was substantially 
constant, though fact varied from 19.6 20.3% the feed. 
This variation not sufficient affect the rate 
The average values each bed temperature are summarized 
Table 

activation energy for the overall oxidation ethylene 
12.9 cal./mole was obtained the usual methods. 


Oxygen Dependence 


The dependence the overall oxidation rate oxygen was 
obtained constant W/F, and constant ethylene concentra- 
tion varying oxygen nitrogen ratio. Data are listed 
Table along ‘with the initial rates calculated directly from 
the conversion data. When log plotted against 
log p., substantially linear were with slopes 
0.23, 0.27 and 0.24 237, 250, and 270°C. respectively. Using 
average value 0.25 for the slope, and combining this result 
with (3), the initial rate overall ethy lene consump- 
tion expressed 


From the data Tables and II, has been evaluated each 
12,900 


Effect Total Pressure 


Data for the pressure range 4.5 atmosphere are listed 
Table Cylinder nitrogen, and oxygen were used place 
air these experiments. Some ‘the data are plotted 
Figure showing that constant total conversion, 
rapidly, then more slowly the total pressure was 
increased. 

Observed rate data finite conversions are lower than 
would predicted Equation This due part the 
inhibiting effect the products the reaction rate. 
has shown that inhibits the rate oxidation ethylene 
ethylene oxide with significant effect the rate oxidation 
ethy lene carbon dioxide and water. reports apparent 
effect either reaction rate when water vapor added with 
inlet found that all three 
products inhibit both main reactions. The order decreasing 
inhibiting strength given ethylene oxide, carbon dioxide 


and water. The inhibiting effects products and the effect 


total pressure were taken into account correlating the data 
with equation the form 


T Peo 
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TABLE 
CATALYTIC OXIDATION ETHYLENE WITH AIR ATM. 


gm. 
3.14 3.81 4.66 6.39 
3.50 4.37 5.37 7.12 
1.41 1.74 2.15 2.81 
*W/F catalyst/mole/hr. total gas flow 
mole consumed/hr. —gm.Ag. 
the selectivity values this group data narrow range 


60, the partial pressure carbon dioxide even the sum 
products may just suitable. 

integrated form this empirical rate equation was used 
for evaluating was obtained substituting Equation (6) 
the flow Equation (7) for piston flow. 


kop !4 

The constant may calculated substitution the experi- 
mental values shown Table into Equation (8). The cal- 
culated values are shown Table These values were 
plotted log 1/T, and show considerable scatter, but 


definite trend with temperature. the method least squares 
2.16 
With the following substitutions for partial pressures, 
TABLE 
and assuming the partial pressure oxygen constant, the atm. 
239 1.68 .086 158 .089 
TABLE 242 2.36 .098 170 
EFFECT VARYING PARTIAL PRESSURE OXYGEN 4.53 115 
OXIDATION ETHYLENE 250 1.0 .093 111 .090 
.063 .088 3.38 4.60 181 175 
270 .187 .750 0.63 5.50 4.60 220 
mole 
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Figure 2—Total conversion C,H, W/F several 
temperatures and constant air:C,H, ratio 15:1. 


Calculated values from Equation (9) are given 
Table These values along with the empirically deter- 
mined values and were used calculate the conversion 
Conversion values Table and the solid lines drawn 
Figures and which were calculated from Equation (8) show 
satisfactory agreement with the observed data. 


Selectivity 

The selectivity the catalyst showed the usual trend noted 
previous Selectivity increased slightly with 
increasing partial pressure oxygen, and decreased with in- 
creasing catalyst temperature. Data illustrating the trend with 
temperature have been plotted Figure Values selectivity 
were taken from this best fit line calculating the conversion 
curves shown Figures and 


Discussion Results 

The power law form the rate equation was chosen because 
its simplicity. Equations proposed from mechanistic consid- 
erations did not fit these data too well. For example, Kurilenko’s 
equation predicted dependence oxygen concentra- 
tion whereas these data showed some dependence. Orzechowski 
and MacCormack’s Equation (3) which fit their data (at 
atm) similar catalyst was not fully applicable this group 
data. Their equation for initial rate can written the form 


Initial rate data given Table from experiments atmospheric 
pressure fitted Equation (10) and thus were agreement with 


200 220 


TOTAL PRESSURE (ATM.) 


Figure 3—Total conversion C,H, total pressure 
constant W/F 1.64; air: ratio 15:1. 


the earlier This equation also predicted satisfactorily 
the functional dependence initial rate total pressure, but 
the constant which can obtained from the intercept plot 
1/p, was from 1/3 1/18, depending the temperature, 
the value derived from the data atmospheric pressure. 
Thus constants and the higher pressures were not 
consistent with the constants obtained equal one; and 
appear functions total pressure well temperature. 

Power law forms the rate equation have been applied 
this reaction. Cambron and found that the overall 
agreement with the data reported here. correlated 
his data with the equation 


Thus, exponents and vary magnitude and seem 
depend catalyst preparation and conditions under which the 
catalyst has been tested. 

The activation energy obtained from the overall oxidation 
ethylene these tests 12.9 cal./mole. This value 
somewhat lower than given Kurilenko who reported 
15.2 and 19.8 cal. for reactions (1) and (2) respectively. 
The reason for the lower activation energy not known, 
though might possibly due the onset diffusion effects. 
rough estimate the eddy diffusivity coefficient within the 


ATM. 


260 


TEMP. 


Figure 4—Selectivity temperature. 
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fluidized bed using Gilliland and Mason’s data‘® indicated that 
significant radial concentration gradient should appear for 
the reaction rates observed. estimate the back mixing 
effect using the same data, also showed that back mixing should 
have only minor effect the conversion levels. The diffusion 
effects within the catalyst film itself are believed small 
but could not satisfactorily evaluated because effective 
diffusivity data for the catalyst film were not available. 

comparison observed catalyst activity with data derived 
from catalyst coated tubes similar those used 
showed that under comparable operating conditions, the activity 
was substantially the same, but that selectivity may slightly 
lower. Selectivity significantly lower when compared with 
Hayes’ data where conditions were not identical. Hayes used 
lower flow rates, and carefully purified feed gases. 

These experiments have shown that the immersion fixed 
strips inert fluidized bed feasible method obtaining 
close temperature control without appreciable loss catalyst 
performance. 


replacement when becomes necessary. obvious limitation, 
which common fluidized bed systems, the limitation 
maximum flow rate, which thus restricts the space velocity 
relatively low values. This maximum flow rate set the 
flow which severe slugging develops. The scale factors 
which affect this limit have not yet been determined. 


Summary 
Small scale experiments have shown that close temperature 
control without reduction catalyst activity can obtained 
immersing catalyst coated strips inert fluidized bed. 
Rate data obtained from the air oxidation ethylene using 
immersed Cambron silver catalyst coated strips glass bead 
bed were correlated the empirical equation 


po''* 


Constants have been evaluated for variables the range 
230-270°C., p-0-50 p.s.i.g.; ratio 35/1 15/1; 
concentration —56%. The constants and for this 
range data are: 


2.80 


logioks = + 3.822 = 


2.16 


Moreover, the method permits easy catalyst 
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Nomenclature 


oxygen ethylene ratio 


total gas flow mole/hr. 
mole 


empirical rate constants 


mole 

empirical exponents 1/2 and 1/4 resp. 

Peo partial pressure ethylene, oxygen, ethylene oxide, 
resp., atm. inlet gas stream. 

total pressure atm. 

initial rate, and rate consumed mole/ 
hr. gm. catalyst 

gas constant 

selectivity 

temperature 

weight silver catalyst gm. 

fraction inlet ethylene converted 

mole fraction ethylene, oxygen feed 
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Evaporation Rates Spray Drying’ 


Heat and mass transfer coefficients have been 
determined during the evaporation water sprays 
produced pneumatic atomizing nozzles 
vertical, cocurrent spray dryer, 8-in. diameter 
and ft. high, under various operating conditions. 
The progressive evaporation the spray down the 
chamber was followed taking samples the drying 
air and the water droplets various distances 
from the nozzle, and new volumetric method for 
the accurate determination air humidities was 
developed. 


For droplet spray ranging mean diameter from 
11.5 38.5 microns, the heat and mass transfer 
coefficients were found essentially the same 
for single, stationary droplets evaporating still air, 
the similarity being attributed the almost equal 
diffusivity the drop and the air stream, resulting 
negligible relative velocity. Certain important 
aspects spray dryer design are discussed, and 
method outlined for the calculation the residence 
time the evaporation zone. 


drying chamber conventional spray dryer may 
considered consisting three distinct zones, namely the 
nozzle zone, which the newly-atomized drops rapidly deceler- 
ate their final settling velocity, the evaporation zone, 
which the mechanism water vapor evolution similar 
that from free liquid surface, and finally the drying zone, 
which internal diffusion the drying particles becomes the 
governing factor. The present investigation was undertaken 
determine the rate heat and mass transfer the drops the 
evaporation zone under conditions likely encountered 
spray dryers. 


Evaporation from Single Particles 

1910, showed experimentally that the rate 
evaporation from stationary drop gas was proportional 
its diameter. His results were analyzed who 
presented the following theoretical equation: 


which the factor accounted for the thickness the film 
received October 1959; accepted March 30, 1960. 


2Pulp and Paper Research Institute of Canada and Department of Chem- 
ical Engineering, McGill University, Montreal, Que. 
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Les auteurs ont déterminé les coefficients 
transport chaleur masse durant 
d’eau jets pulvérisés produits par des atomiseurs 
pneumatiques, dans séchoir jet vertical con- 
courant, pouces diamétre par pieds 
hauteur. Ils ont suivi progressive jet 
long chambre prélevant des échantillons 
d’air gouttelettes des distances variées 
méthode volumétrique pour déterminer avec pré- 
cision 

Les auteurs ont trouvé que pour des gouttelettes 
ayant des diamétres variant entre 11.5 38.5 
microns, les coefficients transport chaleur 
masse sont essentiellement les mémes que pour des 
gouttelettes uniques stationnaires qui 
dans repos. attribue ces comporte- 
ments semblables aux diffusivités trés voisines 
gouttelette courant donnant une vitesse 
relative négligeable. présente certains aspects 
importants réalisation séchoirs jets pulvé- 
risés décrit une méthode pour caleul 
temps résidence dans zone 


surrounding the drop, through which diffusion the vapor 
occurred, and was defined by: 


Equation (1) was verified many workers 
and was modified others into account certain specific 
conditions. particular, made correction for the 
concentration the diffusing vapor infinite distance from 
the drop, and added term take into account 
the relative velocity between the gas and the drop. His equation 
was based approximate solution the Navier-Stokes 
equations but included experimental coefficients, and can 
written follows: 


Several other forms this equation were suggested, and 
1952 Ranz and following extensive study, pre- 
sented the following four equations for heat and mass transfer 
drop under natural forced convection: 


high rates evaporation which may occur with elevated 
drying air temperatures when the liquid has low latent 
heat vaporization the diffusing vapor leaving the drop 
surface changes the temperature gradient the gas film surround- 
ing the drop, and the apparent the Nusselt Number 
given Equations (4) and (6) are too high. get the actual 
value Nu, use correction gas been recommended 
several 


Evaporation from Clouds Particles 


Although the same theoretical considerations apply clouds 
drops for single drops, the problem complicated the 
fact that, spray dryer, the drops are dispersed turbulent 
gas stream. and studied theoretically the forces 
acting single particles suspended turbulent gas and their 

results indicated that the eddy diffusivities the particles and 

the gas are almost equal for small particles and low 
intensities turbulence. This means that the drops should 
evaporate rate corresponding zero relative velocity. 

studied experimentally the diffusion atomized 
sprays turbulent air streams. The work was carried out 
long duct with air velocities ranging 
from feet per second, and Kesler found that the diffusion 
water and alcohol drops ranging from 
microns diameter were approximately equal that the air. 
also found that the rate evaporation 25-micron alcohol 
drops was equal the predicted rate for stagnant conditions 
and therefore concluded that the relative velocity between the 
drops and the air stream was essentially zero. 

few studies have been reported the rate evaporation 
under decelerating accelerating conditions the nozzle zone 


Figure 1—Experimental spray dryer. 
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spray and particularly the field fuel atomiza- 
Most the data give fair agreement with Equation 
(4), but the high relative velocities between the drop and the 
air tend mask any effects due air turbulence and the results 
cannot readily applied the evaporation zone spray dryer, 


Experimental 


conventional spray dryers the particle trajectories and 
gas flow patterns are invariably quite complex, and therefore 
impossible follow the progressive evaporation the droplets. 
Following the approach initiated special vertical 

cocurrent spray dryer was therefore designed, permitting the 
simultaneous determination the air temperature and humidity, 
well the droplet size distribution evaporation proceeded, 
thus enabling the calculation heat and mass transfer coefficients 
along the drying chamber. 


Equipment 

The experimental spray dryer consisted the same basic 
units industrial installation, except that device for 
product collection was provided. 

The drying chamber was 14-ft. long, 8-in. diameter, and 
was insulated with 2-in. layer vermiculite. Sampling ports 
were provided along the whole length the chamber, shown 
Figure 

Two blowers were used, one each end the unit and both 
with sliding gate, that the air rate and pressure 
the drying chamber could readily controlled. calibrated 
orifice was used meter the drying air and perforated plate 
was installed the top the chamber provide uniform air dis- 
tribution. Three electric heaters with total capacity 6.25 kw, 


Figure 2—Humidity determination apparatus. 
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Figure 3—Photomicrograph water droplets (220x). 


supplied the necessary heat. variac and Aminco bimetallic 
thermoregulator assured fine control the air temperature and 
Speedomax 12-point potentiometer recorder was used 
measure the inlet and outlet air temperatures well the wall 
temperatures various points the drying chamber indicate 
when steady state conditions were established the system. 

internal- -mixing pneumatic nozzle (supplied the Spray- 
Systems Company) was installed the center and the 

the drying chamber. Two different sizes were used, 
(Models No. and 2B), both with spray angle small enough 
insure that the droplets did not impinge walls the 
drying chamber under the conditions studied. The liquid feed 
the nozzle consisted tap water metered through rotameter 
and heated cooled heat exchanger just before entering 
the nozzle. Instruments were provided indicate the atomizing 
air and liquid feed pressures, well the feed temperature. 

All valves, instruments and controls were mounted 
central control panel, portion which can seen the bottom 
Figure 


Methods Measurement 
(a) Air Temperature 


The presence liquid droplets the air stream necessitated 
the use special shield protect the temperature measuring 
device. After considerable experimenting was found that 
semi-circular shield, 0.4-in. diameter and 1.5-in. long, pro- 
tecting mercury thermometer bulb 0.2-in. diameter, was 
most suitable. spraying colored solutions into the chamber 
any droplet impingement the bulb could observed after 
withdrawi ing the thermometer and was found negligible. 
Conduction losses along the stem were reduced installing 
the thermometer glass tube. 


(b) Air Humidity 

Owing the small changes drying air humidity observed 
this investigation, considerable degree accuracy was 
required for the humidity measurements. The majority the 
more conventional methods were found either insufficiently 
accurate unsuitable because large air sample requirement. 
volumetric method determining the humidity the air 
was therefore developed, based the absorption water 
vapor magnesium perchlorate (which does not absorb carbon 
dioxide any other common gases) under constant temperature 


and pressure conditions, and the accurate measurement 


the resulting change the volume the air sample. 

precision Mine Air Gas burette 100-ml. capacity, 
calibrated 0.05-ml. subdivisions the range 94-100 was 
used measure changes volume the air low 0.01-ml. 
after was passed over magnesium perchlorate U-tube. 
From the burette the air could passed through the U-tube 
and into self-sealing mercury reservoir, consisting two 
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interconnected chambers, mounted one above the other. The 
air entering the lower compartment forced the mercury into 
the upper section. When the air was sucked back into the burette, 
the mercury flowed back into the lower chamber. fixed level 
mercury the lower chamber was obtained means 
pilot-light circuit with contacts sealed into the mercury chamber. 
passes were found required for complete removal 
the water vapor. 

constant temperature water bath, and jacket around the 
burette with water circulating from the bath through it, wer 
used maintain rigidly constant temperatures throughout 
system, while compensator tube and manometer ensured 
constant reference pressure. photograph the apparatus 
shown Figure 

The whole method depended the assumption that the 
volume the magnesium perchlorate did not change during 
test that the volume dry air trapped the U-tube was 
the same for the initial and final measurements. This assumption 
was undoubtedly justified view the small amounts water 
which were absorbed, but the equipment was calibrated 
ensure Air known humidity was prepared 
saturating given temperature two absorption columns 
which can seen the right Figure Calibration tests 
indicated that the method was accurate within 0.78% the 
range 0.0118 0.0150 air. 

Air samples were obtained from the chamber employing 
inverted cup, 0.75-in. which was inserted into the drying 
chamber and through which air was sucked into evacuated 
sampling bottle. The rate suction was such that the air 
velocity the cup was well below the average terminal velocity 
the drops, thus achieving the desired separation. aluminum 
lip the periphery the cup prevented the collected spray 
from entering the sampling device. 


Droplet Size 


Since the rate evaporation function the droplet 
diameter, was necessary determine the particle size distribu- 
tion various locations the drying chamber. The microscopic 
method, being direct, was chosen this inv estigation. 


Samples the spray were collected immersion cell 
filled with Varsol, similar that described 
prevent aporation the droplets before they could photo- 
graphed. The immersion cells (0.18-in. 1.D. and 0.25-in. O.D.) 
were made minimize deflection the small particles, 
and piece optical glass, coated with General SC. 
was cemented into the bottom the cells. Rupe 
has shown that, employing the above precautions, the drops 
remain spherical and their diameter can therefore measured 
directly. Furthermore, the low surface tension between Varsol 
and water keeps coalescence between the drops 
Good target efficiencies were obtained with these small cells, 
evidenced the large number small drops (down 
collected. 

The sample obtained from traverse the spray was 
immediately photographed through microscope with 
magnification using single lens reflex camera and high resolu- 
tion film (Adox KB-14, 141 lines per mm.). High contrast 
prints were obtained from the negatives and these were viewed 
and counted under low magnification through micro- 
scope. Several hundred drops were measured from each sample 
and grouped into sizes separated usually 2.38 microns. 
typical photomicrograph water droplets shown Figure 


Procedure 

series runs was made determine the rate evaporation 
the water droplets under various experimental conditions, 
according rigidly standardized operating 

After selecting the proper air temperature and air rate, the 
system was allowed reach steady state conditions. feed 
was introduced the nozzle this stage, but atomizing air 


115 


was passed through the desired pressure. 
traverse the air along the length the drying chamber was 
then made determine the heat losses and correct for the 
presence the cold atomizing air. was invariably found 
that the heat loss was negligible, but that the correction for 
atomizing air was necessary corrected inlet air temperature 
thus obtained was designated 


Water was then introduced the nozzle and once again 
time was allowed for steady state reached. Special pre- 
cautions were taken this stage ensure that the spray did 
not impinge the chamber alls shown observation 
the wall temperatures and physical inspection the walls with 


proper lighting. 


The feed, atomizing air and inlet drying air conditions and 
rates were noted and the temperature and humidity distribution 
along the drying chamber were determined described pre- 
viously. Droplet size distributions were also obtained selected 
distances from the nozzle. was ascertained, means 
unshielded thermocouple probe, that all runs the water drop- 
lets remained the wet bulb temperature the air throughout 
the chamber. 


Caleulations 


The rate evaporation the water droplets was calculated 
directly from either heat balance material balance. For 
the accurate calculation the instantaneous heat transfer co- 
efficients, however, equation had derived differential 
basis because the continuous variation the values the 
variables involved the system. Considering height 
the chamber, the following heat balance can written: 


while the rate heat transfer the drops located the 
differential volume the chamber height given by: 


RUN 
RUN 
RUN 


160 


150 


DISTANCE FROM NOZZLE 


Figure 4—Typical temperature distribution spray dryer. 
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(9) 


the surface area the droplets can obtained from 
the specific area the spray (in water), means 
the expression: 


Equation (10) the terminal velocity the droplets, 
considered negligible comparison with the absolute 
velocity the air, the actual experiments, varied 
from ft. per sec., depending the drop diameter, 
the pounds unev aporated liquid droplets passing per hour 
through the cross-section question, was calculated from the 
water feed rate the nozzle means humidity balance, 
and was obtained from the droplet count the section 
substitution into the following expression: 


From these equations, the instantaneous values could 
readily calculated. These are the actual values from which the 
actual Nusselt Numbers, can obtained. The 
latter can compared directly with the apparent 
dicted Equation (4) since the drying air temperatures 
(maximum 229.7°F.) were low, and the correction for high 
mass transfer rate consequently negligible. 

From similar equations and from their definition, the instant- 
aneous mass transfer coefficients, and the modified Nusselt 
Numbers for mass transfer, Nu’, were calculated. 

noted that the mean film properties such 
py, ky, were estimated the arithmetic average temperature 
between the air and the surface the droplet, the air wet 
bulb temperature. 


FRACTION WATER VAPORIZED 


SECONDS 


Figure 5—Rate water vaporization. 
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Figure 6—Variation heat transfer coefficient with mean 
statistical droplet diameter, 


Results 


total ten runs were made under different conditions 
inlet air temperature (135.2 229.7°F.) air velocity (3.9 
14.8 ft./sec.) and average droplet size (11.5 38.5 The 
initial average droplet size was controlled either changing 
the atomizing air pressure using different nozzle. All 
runs were made under pressure and effort was 
made control the inlet air humidity, which ranged from 0.0054 
0.0123 lb. water dry air. 
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Figure 7—Comparison experimental data with theoretical 
relation 


Some the observed values the experimental variables 
well the calculated results are shown Tables and 2*. 
The calculated results include the amount water passing any 
given section calculated from the heat balance and 
material balance and the slope the curves the 
temperature versus distance from nozzle (dt/dx). 

Three typical longitudinal temperature distributions along 
the spray dryer are shown Figure Figure shows the 
fraction water vaporized function time the spray 


TABLE 
HEAT AND MASS TRANSFER WATER DROPLETS LARGE AVERAGE DIAMETER 
Nozzle No. 22B Atom. air pressure Feed temp. 95°F 
Feed rate 6.42 Drying air rate 309 
Average air velocity 3.9 ft./sec. 95°F. 
OBSERVED CALCULATED 

0.79 168.0 38.5 3.22 270 2.05 1.91 

154.7 36.9 2.33 —31.0 299 2.20 2.01 

1.46 145.5 1.62 23.0 328 

1.79 138.3 0.0255 30.8 1.07 343 2.14 1.94 

1.96 0.97 
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200 


TABLE 


HEAT AND MASS TRANSFER TO WATER DROPLETS OF SMALL AVERAGE DIAMETER 


no. 
Nozzle No. Atom. air pressure p.s.i.g. Feed temp. 85°F 
Feed rate 4.86 Drying air rate 860 
Average air velocity 10.7 ft./sec. 86°F 
OBSERVED CALCULATED 
0.96 154.0 18.3 25.7 560 2.06 1.92 
1.13 150.3 2.71 21.1 571 2.03 1.88 
1.29 147.1 16.9 610 2.09 1.92 
1.46 144.2 15.7 13.8 660 2.09 1.94 
1.79 140.5 13.6 0.82 0.90 9.0 737 2.04 1.86 
1.96 139.2 126 6.3 760 1.96 1.79 
2.54 137.0 11.5 0.12 863 2.02 1.85 
dryer. The very striking effect droplet size the rate Conclusions 


evaporation can seen from these curves: for instance, runs 
No. and were carried out under identical conditions, except 
for the initial droplet size. 

Finally, the heat transfer were plotted Figure 
against the mean Sauter diameter, while Figure shows 
comparison between the same experimental points and straight 
line corresponding theoretical Nusselt Number 


Discussion Results 


The experimental procedure followed this investigation 
was planned permit complete analysis the process 
evaporation from water droplets suspended turbulent air 
stream. ‘The good agreement between the heat and material 
balances and the Tables) indicate the general 
accuracy the results. 

The progressive evaporation the water droplets along the 
drying chamber clearly represented the temperature 
distribution curves, typical examples which are shown 
Figure Starting with gentle slope the immediate vicinity 
the nozzle, the slope slowly increases maximum, all 
curves exhibit point beyond which 
exponentially until their slope nearly zero. detailed study 
the nozzle range has already been presented 
and confirmation his observations the present tempera- 
ture distribution data indicate that the amount water evap- 
orated the first few inches from the nozzle comparatively 
This due the fact that the droplets leave the nozzle 
high initial velocity, and although their rate evaporation 
high (due their high relative velocity) their residence time 
the nozzle zone too permit any large amount 
evaporation. 


the evaporation zone shown Figure the average 
value the Nusselt Number for all experimental points was 
found 2.07 The average value the modified 
Nusselt Number, Nu’, was 1.89 0.06, but this somewhat 
lower value can attributed uncertainties the estimation 
the diffusivity 

*Tables 1 to 10 of this paper have been deposited as Document No. 
6344 with the ADI Auxiliary Publications Project, Photoduplication Ser- 
vice, Library of Congress, Washington 25, D.C. A copy may be obtained 
by citing the Document No. and by remitting $1.25 for photoprints, or 
$1.25 for 35 mm. microfilm. Advance payment is required. Make cheques 
or money orders payable to: Chief, Photodupliaction Service, Library of 
Congress. 
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The results obtained indicate conclusively that the rate 
evaporation spray drying can safely calculated assuming 
that the individual droplets evaporate rate corresponding 
stagnant conditions, i.e. assuming the Nusselt Number 
equal two, predicted Equation (4) for This 
agreement with the results who also showed 
that, for drops the same diameter range the present 
study, the effective relative velocity between the droplets and 
the air stream was pratically zero, spite the air turbulence. 
furthermore apparent that the presence the droplets 
the form cloud has influence the rate evaporation, 
least the voidage ratio commonly encountered spray 
drying. 

Because accounts for the bulk the evaporation which 
occurs spray dryer, the evaporation zone considerable 
importance the development theory spray dryer 
design. calculate the residence time needed vaporize the 
largest drops produced the nozzle which will determine the 
length the evaporation zone for given air flow itis 
necessary estimate the progressive decrease air tempera- 
ture and humidity evaporation progresses, since the latter 
dictate the instantaneous rates heat and mass transfer the 
largest drops. The calculations can readily carried out 
iterative method over small increments time, the drop 
size distribution known. proposed the 
method consists dividing the spray into small number 
drop size groups (from six ten, depending the accuracy 
required) and performing the calculations each size group 
over small increment time, using the d,, for that group 

The accuracy the method demonstrated Figure 
which compares the experimentally observed size distribution 
run No, distance 5.2-ft. from the nozzle, with that 
calculated from the experimental drop distribution measured 
3.19-ft. from the nozzle. agreement between the two 
distributions excellent. 

Once the evaporation occurring each size group over 
small increment time has been determined, easy 
matter calculate the resulting changes air temperature and 
humidity for the interval. turn permits calculate the 
cor responding evaporation from the largest drops produced 
the atomizing nozzle. 

‘The method approach which has been outlined requires 
knowledge the initial drop size distribution. the latter can 
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Figure 8—Calculated and experimental droplet size 


distribution. 


expressed the form simple mathematical probability 
tion and the vapor pressure the solvent can assumed 
remain unaffected the presence solute, the analytical 
method Probert can used, and tedious iterative solution 
necessary. (This method, however, ignores the tempera- 
ture drop the air evaporation proceeds). number 
such distribution functions have been proposed for certain atom- 
izing general, however, will necessary 
determine the distribution (and the size the largest drop) 
experimentally for given nozzle and for specific feed liquid. 
The small amount evaporation occurring the nozzle zone 
can either neglected (which will lead conservative 
estimate the evaporation zone proper) or, preferably, the 
used. completely analytical solution the complex mass 
and heat transfer phenomena occurring the nozzle range will 
only possible when accurate data for the coefficient drag 
that region become available. Such study currently 
under way. 


The present investigation has been concerned with the evap- 
oration clouds water alone. The information obtained can 
easily extended more complex cases, some which have 
already been studied, such the effect the presence solute 
the feed solution and the behavior droplets evaporating 
Many more considerations still require 
investigation, particularly those associated with the turbulence 
aspects the system, but can now predicted that solution 
the problem spray dryer design sound theoretical 
basis may soon possible. 
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Nomenclature 


transfer area, 


Cy, capacity gas film surrounding drop, (B.t.u.)/ 

heat capacity water vapor, (B.t.u.)/(Ib.)(°F.). 

humid heat capacity, 

diffusivity water vapor air, ft.?/hr. 

individual drop diameter size distribution count, 
microns. 

d,, mean drop Sauter diameter, ft. 

acceleration due gravity, ft./hr.? 

heat transfer coefficient, 

air humidity, dry air. 

mean thermal conductivity the gas film surrounding 
drop, 

rate liquid spray, calculated from heat 
balance; calculated from material balance. 

M,, average molecular weight gases between the drop 
surface and the bulk air, 

number droplets size distribution count. 

water vapor pressure bulk air, 

mean pressure air film, p.s.f.a.; defined by: 

equilibrium vapor pressure the drop surface (corre- 
sponding p.s.f.a. 

radius droplet, ft. 

outer radius film around droplet, ft. 

specific surface area drops, 

temperature air, °F. 

wet bulb temperature air, °F. 

velocity drying air relative drop, ft./hr. 

absolute velocity drying air, ft./hr. 

terminal velocity drops, ft./hr. 

distance from nozzle, ft. 

mean absolute viscosity gas film surrounding drop, 

atmospheric pressure, p.s.f.a. (in definition p.). 

time, hr. 

References 

(1) Morse, W., Proc. Am. Acad. Arts and Sci., 45: 363 
(1910). 


(2) Langmuir, I., Phys. Rev., 12: 368 (1918). 

(3) Fuchs, N., Physik. Sowjetunion, 224 (1934). 

(4) Froessling, N., Gerlands Beitr. Geophys., 52: 170 (1938). 

(5) Ranz, E., and Marshall, R., Jr., Chem. Eng. Progr., 
48: 141, 173 (1952). 

Symp. Combustion”, Cambridge, Mass. (September 
1952). 

(7) Marshall, R., Chem. Canada (1955). 

(8) Spalding, B., Proc. Roy. Soc. (London), 
(1954). 

(9) Ranz, E., Trans. Am. Soc. Mech. Eng., 78: 909 
(1956). 

(10) Sleicker, A., Jr., and Churchill, W., Ind. Eng. Chem. 
48: 1819 (1956). 


119 


/ 


(11) Liu, Vi-Cheng, Dept. U.S. Air Force, Project 2160 
(1955). 

(12) Soo, Chem. Eng. (1956). 

(13) Kesler, H., Sc. Thesis, Mass. Inst. Technol., (1952). 

(14) Manning, P., and Gauvin, H., (In 
press). 

(15) Ingebo, D., Nat. Adv. Comm. Aero. Tech. Note 3762 
(1956). 

(16) Fledderman, G., and Hanson, R., Eng. Res. Inst., 
University Michigan, Report No. 667 (1951). 

(17) York, L., and Stubbs, E., Trans. Am. Soc. Mech. 
74: 1157 (1952). 

(18) Pinder, L., Eng. Thesis, McGill University, Mon- 
treal, Quebec (1952). 


(19) Rupe, H., “Third Symposium Combustion 
and Explosion Phenomena”, The Williams and 
Co., Baltimore (1949). 

(20) Marshall, Jr., Trans. Am. Soc. Mech. Engrs., 

(21) Probert, P., Phil. Mag., 37: (1946). 

(22) Marshall, R., “Atomization and Spray Drying’, 
Chemical Engineering Progress Series, No. Vol. 

1960). 

(24) Gauvin, H., Knelman, H., and Pinder, 
presented the Symposium Spray Mechanisms, 
nual Meeting. Cincinnati, Ohio (Dec. 1958), 


The Canadian Journal Chemical Engineering, August, 1960 


LETTE 


Vap 
C.J. 


Editor. 
The 
enables 
li uid 
solving 
two 
other 
also 
and 
work). 
there 
abou 
are 
are 
are 
though 
33. Its 
One 
claim 
log 
even 
Coull 
(1952)) 
and the 
their 
The 


Editor 

like 

integra 
intende 
tion 
equatio 
their 
restrict 
the 
seems 
data 
the 
method 


The 


| 


THE EDITOR 


Note Thermodynamic Consistency Ternary 
Vapor-Liquid Equilibrium Data, and Lu, 
37, 117-120, 1959. 


London, England. 
Editor, 


The authors have discussed new way applying estab- 
lished relationships experimental data, and claim that 
enables them select inconsistent values from set vapor- 
liquid equilibrium data ternary system. This claim can dis- 
puted two regards. Firstly, for the three methods suggested for 
solving their integral Equation (5), they have only provided 
two examples calculations. This itself highly selective 
and particularly unwarranted comparison not made 
other methods previously used for evaluation the same 

the second place, within their calculation sets there 
also unwarranted selection data. For example, Table 
and Figure ten experimental points from their reference are 
considered (corresponding runs 29, the original 
work). The authors then test these points and show that 
there “overall deviation thermodynamic consistency” 
about 38%, which reduced points 29, and 
are omitted. They therefore claim that these latter points 
are inconsistent. However, the original points plus run 
are examined, the deviation drops from 38% below 14% even 
though includes the so-called inconsistent points 29, and 
33. seems clear then that the method advocated the authors 
neither sensitive nor critical enough meet the claims made. 

One other minor criticism may also made. The authors 
claim that the testing isobaric data, the assumption 
Broughton and Brearley (Ind. Eng. Chem. 47, 838, 1955) that 
even for binary The basis for their claim and 
Coull (Chem. Eng. Progress Symposium Series 48, No. 
(1952)). This reference discusses only systems this context, 
and then only from the standpoint White’s equations (Trans. 
41, 539, 1945). Broughton and Brearley validate 
their claim reasonably well for binary and ternary systems. 

The authors approach may well promising one but 
extensive testing will required. 


Shemilt, 


Department Chemical Engineering, 
University College, London. 


Ottawa, Ont. 


Editor, C.J.Ch.E.: 


wish thank Dr. Shemilt for his kind interest our 
publication. reply his criticism and comments, would 
like offer the following: 

our paper, the proposed graphical and numerical 
integration methods are developed theoretically. They are 
intended illustrate the usefulness the Gibbs-Duhem equa- 
tion its exact form. This not the case where empirical 
equations are proposed, using data fitting the only means for 
their justification. the development the proposed methods, 
involved other methods are avoided, indicated 
the manuscript. provide too many examples calculation 


seems wasteful and unnecessary. The idea extensive testing 


data for supporting the Gibbs-Duhem equation untenable and 
the validity exact method should not judged empirical 
methods simply because they are available. 

might mentioned that our proposed methods have 
already been employed the literature for testing ternary 
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Figure 1—Arrangement experimental points for thermo- 
dynamic consistency test. 


vapor-liquid equilibrium data (e.g. Free and 
Hutchison, Chem. and Eng. Data 306 (1959)). 

wish point out that the calculation made the 
criticism does not seem correct. If, from the Reference 
(8) the manuscript (John Griswold and Wong, Chem. 
Eng. Progress Symposium Ser. 48, No. (1952)), the ten 
points corresponding runs 29, 39, plus run for the 
system acetone-methanol-water 100°C. are examined, the 
logical sequence these points should arranged indicated 
Figure seen Table that the overall deviation 
the thermodynamic consistency these eleven points 32.3%. 
the calculation repeated, excluding points 29, and 33, 
the deviation reduced 2.4%. The same conclusion 
reached that presented the manuscript. unfortunate 
that the criticism misleading. have not included point 
our paper because little too far away from set 
closely spaced points, required the numerical method. 
should mentioned, however, that not claim that there 
should possibility mutual cancellation positive and 
negative errors. Should the percentage deviation series 
points changed very considerably including excluding 
single point, would obvious that this point thermo- 
dynamically inconsistent. 

the article, mentioned that the validity the 
independent temperature, questionable. still believe 
this and share the point view and Bennett (Ind. 
Chem. 48, 679-80, (1956)) regarding the Broughton and 
Brearley method (Ind. Eng. Chem. 47, 838, (1955)). 


where the partial molal heat solution component 
The assumption that constant composition independ- 
ent temperature, employed Broughton and Brearley, 
indicates constant partial molal heat solution, which obvi- 
ously not necessarily true. Furthermore, requires ideal 
entropy mixing for all solutions. This shown the 
following expressions: 


OGE 


ns 

ky 
Q39 
\ 
- 
037 4 
33 
30 32 70 
= 


TABLE 
Acetone-Methanol-Water 100°C. 


Griswold and Wong, Chem. Progress mposium Ser. 48, No. (1952)) 
log log log 
0.119 0.5035 —0.01297 0.119 0.0170 0.762 0.0785 
0.089 0.5826 —0.01064 0.131 —0.0110 —0.00665 0.780 0.0828 
0.066 0.6230 —0.01155 0.119 0.0678 —0.02111 0.815 0.0588 .03994 
0.024 0.7576 —0.00006 0.107 0.1664 —0.00728 0.869 +0.01790 
0.043 0.5110 0.225 +0.01140 0.732 0.0678 —0.02938 
0.069 0.5136 +0.00128 0.198 0.733 0.0795 —0.02082 
0.100 0.4925 0.193 0.0318 0.707 0.0962 —0.01179 
Repeat the above calculation without run nos. 29, and 
the literature and reported that they vary linearly with tempera- 
ture for many systems over moderate temperature range. 


where and are the excess free energy and entropy 
mixing, respectively. Such solutions having ideal entropy 
mixing are classified “regular” solutions. invalidity 
the assumption that all solutions behave regular solutions 
beyond argument. For detailed discussion, see for example, 
“Chemical Thermodynamics’ Prigogine and Defay, trans- 
lated Everett, Chapters XXIV and XXV, Longmans Green 
and Co. London, (1954). 


1956, Smith and Bennett (Ind. Eng. Chem. 48, 680, 
Column lines 32-4, (1956)) their review article thermo- 
dynamics pointed out that assumption that log 
constant constant composition only roughly true few 
cases and often Recently, (Can. Chem. 
37, 193, (1959)) correlated some heats mixing data available 
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other words, the quantity log constant composition should 
involve three terms, the form (b/T) log 
addition the and Coull (Chem. Eng. Progress 
Symposium Series 48, No. (1952)), one may cite many 
other evidences show that the assumption, that log 
constant composition independent temperature, indeed 
questionable. 

James Li, 

Monroeville, Pa., 


and 


Benjamin C.-Y. Lu, 
Chemical Engineering Department, 
University Ottawa. 
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INDUSTRIAL SECTION 


Effect Diesel Locomotive Operation 
Atmospheric Conditions 
Railway Tunnel’ 


extensive study was made over seven day 
period determine would feasible, from 
air pollution view point, operate diesel loco- 
motive powered trains through the St. Clair Tunnel. 
The test involved the passage approximately 200 
diesel powered freight and passenger trains through 
the tunnel; 3,000 air samples were taken and meas- 
urements were made air flow and smoke density. 
Concentrations carbon monoxide, nitrogen dioxide, 
total oxides nitrogen, aldehydes, sulphur dioxide, 
carbon dioxide and oxygen were determined. 

The sampling methods, sampling equipment, ana- 
lytical procedures and test equipment used for this 
particular investigation should applicable 
variety air pollution problems met the chemical 
industry. 


arrival the era the diesel locomotive the Canadian 

National Railways made confirmation necessary that diesel 
exhaust fumes low concentrations were not detrimental 
atmospheric conditions various installations. 

One the early developments the change diesel power 
required this confirmation the case the replacement 
electrified dieselized trains the Sarnia Tunnel, which forms 
major international traffic link betweenSa rnia, Ontario and 
Port Huron, Michigan. 

This paper, therefore, presents the engineering and chemical 
results tests carried out this tunnel under fully dieselized 
operation October 1957 and October 1958 and deals with the 
techniques approaching and analyzing the factors involved 
the investigation. The toxicological aspects the work have 
been presented paper the Industrial Health Conference 
the American Industrial Hygiene 

Similar tests have been carried out other but 
far known, the present work represents one the most 
intensive studies railw tunnel environment under condi- 
tions fully dieselized operation which has been carried out 
date. For the purposes the tests, was necessary 
suspend all electrically hauled traffic and replace diesel 
powered trains for period seven day Chemical, engineer- 
ing and medical staffs were involved the inv estigation. 

test data obtained confirmed that exhaust fumes result- 
ing from full dieselization the through the tunnel were 
not significantly detrimental atmospheric conditions the tun- 
nel. mean concentrations constituents exhaust gases 
general, below the threshold limits established the 


\Manuscript received September 10, 1959; accepted February 12, 1960. 


Research and Development, Canadian National Railways, 
Montreal, Que. 


National Health and Welfare, Ottawa, Ont. 
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American Conference Governmental Industrial Hygienists. 
There were abnormal health conditions observed 
ticipants during hr. exposure the tunnel atmosphere. 

order reduce concentrations smoke and total oxides 
nitrogen well below acceptable levels and provide 
margin safety the event unusual operation conditions, 
was recommended that system forced ventilation 
installed before full dieselization was introduced standard 
operating practice replace electrical operation trains. 


Diesel Engine Exhaust 


Diesel exhaust fumes contain certain gases resulting from 
the complete partial combustion diesel fuel oil and also 
formed during combustion with the excess air supplied engine 
The more harmful components the exhaust 
gases are carbon monoxide, carbon dioxide, oxides nitrogen, 
oxides sulphur, formaldehyde and The concentra- 
tions which these constituents occur diesel exhaust largely 
depend the fuel/air ratio which turn varies with the load 
the 

The trials the St. Clair Tunnel were conducted using 
1750 Class GR-17 road switchers General Motors manufac- 
ture. these engines the air blower, providing air for combus- 
tion, driven directly off the drive shaft and therefore air 
intake proportional engine speed. Figure shows the 
volumes combustion air, exhaust gases and excess air 
various engine speeds for units this type. 

given speed, the weight air intake per unit time 
approximately constant, whereas the weight fuel injected 
increases the load the engine increases. Under normal 
working conditions the excess air used may much greater 
than the stoichiometric requirement. 

study carried out 1946 Berger and 
reported findings diesel exhaust gas 
locomotive investigated developed 1570 h.p., maximum speed, 
from each its four units. Apart from these results, there has 
been little published information the exhaust high horse- 
power locomotive engines such those used throughout these 
tests. The Canadian National Railw ays Research aboratories, 
however are conducting inv estigations this field, the 
which will reported future date. Preliminary results 
indicate that the range constituent gas concentrations the 
exhaust gas new reconditioned engines 1750 h.p., 
shown Table 


The St. Clair Tunnel 


The St. Clair Tunnel has length ft. and its 
midpoint lies ft. beneath the bed the St. Clair River. 
iron shell ft. diameter and having cross- 
sectional area 312 sq. ft. and volume 1,875,000 cu. fr. 
(Figure has single track and the slope the approach 
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TABLE 


CONCENTRATIONS SOME EXHAUST CONSTITUENTS 
FROM NEWLY RECONDITIONED 1,750 H.P. LOCOMOTIVES 


Concentration p.p.m. 


Exhaust 
Half Load Full Load 

| | | 


and outbound grades 2%. Because the unique elevations 
the tunnel and the cooling effect the river and the mass 
soil above the tunnel, was anticipated that there would 
little, any, natural ventilation. Nevertheless, view the 
operational train speed m.p.h. and the fact that the frontal 
area diesel car unit approximately one-half that the 
cross-sectional area the tunnel, piston action the train 
could expected draw certain quantity diluting air into 
the tunnel with beneficial 


Methods Sampling and Analysis 


The shape, dimensions and the absence forced ventilation 
the tunnel suggested the possibility variations exhaust 
gas concentrations different points along the length the 
tunnel. addition, trains enter the tunnel down grade 
under their own momentum with engines idling and coast 
about the mid-point; here power gradually applied until the 
locomotive engines are full power the outbound grade. 
order contend with this situation, was necessary collect 
air samples representing not only average conditions, but also 
fluctuating, peak and residual conditions. 

Sample collecting stations were therefore chosen points 
600 ft. inside the tunnel portals and the mid-point the 
tunnel insure that the samples collected would represent- 
ative tunnel conditions. diagram showing the locations 
the various test stations shown Figure 

Personnel engaged the work air sampling were em- 
ployed hr. 3-shift basis. all, 218 trains passed 


through the tunnel from both easterly and westerly directions 
during the test. Trains were powered with from one four 
diesel units and speeds were varied between and m.p.h. 
Special idling tests inside the tunnel were conducted using two 
unit diesel locomotives. 


THROTTLE POSITION 
275 275 35 435 515 - 595 
Pp 


Figure 1—Intake air and exhaust. 
1750 h.p.-G.M. engine. 


AIR SAMPLING 
STATIONS 
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air samples were collected automatic abso 

tion samplers which measured volume air controlled 

calibrated orifice passed through the absorption solution for 

minutes. Average samples were collected for nitrogen 

dioxide, carbon dioxide, sulphur dioxide and aldehydes. 
air samples were collected breaking the 

seals calibrated one-litre glass tubes that had been evacuated 

samples were collected for carbon monoxide, nitrogen dioxide, 

total oxides nitrogen, aldehydes and oxygen. The glass inlet 

seals were snapped off the following specified times: 

Summit Train leaving yard 

Zero time Train passing sampling station 


Continuous sampling for smoke and particulate matter was 
carried out Hemeon Samplers and High Volume 
Samplers 

Special sampling for blood tests personnel the three 
gas sampling stations the tunnel before and after each hr, 
exposure was conducted Medical Staff the Department 
National Health Welfare and the Canadian National Railway 
Company. These blood samples were subsequently analyzed 
for carboxyhaemoglobin content. 

Air flow measurements, temperature, relative humidity and 
barometric pressure observations were made during and after 
the passage each test train. 

Methods analysis were chosen, taking into account factors 
such possible comparatively high very low occurrences 
gas concentrations, well the facilities which would 
available the temporary laboratory Sarnia. Consideration 
was also given the effects storage and transportation 
samples which might have sent some 400 miles the 
company’s main laboratories, for final analysis. 

Nitrogen dioxide and aldehydes were determined soon 
possible (approximately hrs.) after collection the samples. 

Nitrogen dioxide concentrations were determined the 
Saltzman colorimetric procedure. This method depends the 
reaction nitrogen dioxide, sulphanilic acid and 
ethylenediamine acetic Sensitivity the method 
about 0.05 p.p.m. 

Aldehydes were absorbed aqueous solution phenyl- 
hydrazine hydrochloride followed treatment with potassium 
ferricyanide alkaline The sensitivity this 
procedure approximately 0.2 p.p.m. 

Sulphur dioxide was absorbed dilute solution sulphuric 
acid and hydrogen Specific conductance the 
resultant sulphuric acid was measured conductivity bridge. 
The particular instrument used was supplied the Industrial 


EASTERN 
ENO 


SARNIA 
PORTAL 
SUMMIT 


PORT HURON 


n 
| SARNIA 


ON. 


‘ATION 9 


STA’ 
STATIONG 
STATIONS 
STATIO 
STATION 3 


STATION 
(FIELO 


Figure section St. Clair tunnel and 


location air sampling stations. 


cell 
ity 
the 
Carbo: 
barium 
amount 
this meth« 
Total 
were 
repared 
This 
passing 
ing 
palladium 
standards 
0.074% 
0.2 
continuou 
The worl 
urement 
carbon 
this ins 
Carbo 
means 
Roughlan 
measured 
smoke 
tance. 
and the 
linear fee 
filtrat 
glas web 
determin: 
pressed 
Air 


velomete 


Discussi 


The 
and irrit: 
hr. 
The thre 
most 
toxicity 

The 


concentr: 


WESTERN 
END 
6032' 
5000 
z 
lad 
12 
«2000 
1000 


Company, and was designated Model RC-1. The 
cell the conductivity cell was 0.1 mhos. sensitiv- 
ity 0.05 was obtained these determinations 
the method employed. 

dioxide was absorbed barium hydroxide containing 
barium The concentration was estimated from the 
amount neutralized barium hydroxide. The sensitivity 
this method approximately p.p.m. 

Total oxides nitrogen the time the investigation 
were determined the modified phenoldisulphonic acid method 
prepared the Research Department the Union Oil Company 
California tentative method A.S.T.M. Committee 
‘Total oxides nitrogen were calculated 
This method has sensitivity approximately p.p.m. 

Carbon monoxide instantaneous samples was determined 
passing the air sample through the NBS colorimetric indicat- 
ing gel contains ammonium molybdate and 
palladium sulphate. The color produced was compared with 
standards prepared similar manner from air containing 
0.074% CO. This method has sensitivity approximately 
0.2 p.p.m. Carbon monoxide concentrations were also recorded 
continuously the direct indicating 
The working principle this instrument based the meas- 
urement the heat oxidation reaction carbon monoxide 
carbon dioxide over the catalyst The sensitivity 
this instrument approximately 0.002%. 

Carbon monoxide the blood personnel was determined 
means the microgasometric technique Scholander and 
Some samples were analyzed the infra-red 
absorption 

Smoke concentrations were determined filtration 
measured volume air through paper tape (Hemeon) 
smoke samplers and evaluation the spot optical transmit- 
tance. Smoke was sampled continuously hourly cycle 
and the concentrations were calculated COH units per 1,000 
linear feet air filtered. Particulate matter was determined 
filtration over hr. periods through 1,106 flash-fired Fiber- 
glas web high volume samplers, followed gravimetric 
determination the Fiberglas deposits. The results are ex- 
pressed milligrams per cubic meter air. 

Air flow measurements were obtained means calibrated 
velometers and annemometers. 


Discussion Test Results 


The maximum allowable concentrations each toxic gas 
and irritant the atmosphere for continuous exposure during 
hr. period have been established Committee the 
American Conference Governmental Industrial Hygienists. 
The threshold values are revised each year the light the 
most up-to-date information and experimentally verified data 
toxicity 

The 1958 revision shows the following maximum allowable 
concentrations (M.A.C.): 


Carbon dioxide 5,000 p.p.m. 
Carbon monoxide 100 p.p.m. 
Total oxides nitrogen p.p.m. 
Nitrogen dioxide p.p.m. 


CROSS SECTIONAL 
AREA 312 SQ. FT. 


FREE AREA 


AVERAGE 162 SO. FT 


CAR FRONTAL 
AREA 
150 SQ. FT. 


Figure 3—Cross section St. Clair tunnel. 


Sulphur dioxide p.p.m. 
Formaldehyde p.p.m. 
Typical concentrations the gases found the tunnel 
atmosphere during the passage diesel trains are presented: 


(a) Carbon Monoxide 


Concentrations carbon monoxide found the tunnel were 
well within safe limits. With only exceptions, one which 
can noted Table 150 p.p.m., concentrations were 
well below the M.A.C. 100 p.p.m. for hr. exposure period. 
Table gives summation results the gas sampling station 
the middle the tunnel. 

Carboxyhaemoglobin contents the blood test personnel 
after hrs. exposure the tunnel environment, are shown 
Table The pre-exposure percentages roughly parallel individ- 
ual smoking Post-exposure figures indicate slight 
trend toward higher HbCO saturation percentages approx- 
imately half the subjects. This substantial agreement 
with the fluctuating, but generally low, carbon monoxide 
concentrations found the tunnel atmosphere. (Table 


(b) Aldehydes 


indicated Table showing results for aldehydes 
the sampling station the centre the tunnel, concentrations 
aldehydes were generally below the M.A.C. 5.0 p.p.m. 


(c) Oxides Nitrogen 


Typical values for oxides nitrogen are shown Table 
the gas sampling station nearest the Port Huron portal. 

Average values for all samples taken for nitrogen dioxide 
were well below the M.A.C. p.p.m. 

Total oxides nitrogen expressed nitrogen dioxide found 
approximately 190 samples collected during the passage 
test trains showed concentrations considerably excess 


This value, although not recognized 


TABLE 
MONOXIDE CONCENTRATIONS P.P.M. St. CLAIR TUNNEL STATION (MIDDLE) 


Date Time 
amples 


Highest values exceedin 

Mean Concentration M.A.C. 100 

3.8 

2.7 

2.8 

2.4 
13.6 150 
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TABLE 


CARBON MONOXIDE PER CENT SATURATION 
CARBOXYHAEMOGLOBIN BLOOD PERSONNEL EXPOSED 
TUNNEL ATMOSPHERE For 


Before After 
4.23 4.0 Nil 


nevertheless serves guide possible threshold value, 
since has been established that values below p.p.m. are 
not significant hazard. The average all samples taken was 
approximately p.p.m. with mean values the sampling 
stations varying from 9.1 38.4 p.p.m. and the peak instan- 
taneous concentration reaching level p.p.m. 

This substantial difference between concentrations 
trogen dioxide and total oxides nitrogen results from the 
presence nitric diesel exhaust gas. the high 
combustion temperatures existing diesel engines, the reaction 
2NO takes place. This reaction has been men- 
tioned Davis and Holtz. Temperature and oxygen concentra- 
tions result oxides nitrogen varying from mixtures 
and 

Apart from this, number reactions occur where mixtures 
nitrogen dioxide, nitric oxide, oxygen and water vapor are 
present the diesel These result the conversion 
nitric oxide nitrogen dioxide nitrous acid: 


TABLE 
ALDEHYDES FORMALDEHYDES P.P.M. St. TUNNEL STATION (MIDDLE) 


the equilibrium condition, where the concentrations 
nitric oxide and nitrogen dioxide are equal and about water 
vapor present, the product (nitrous acid) would about 
19% the nitric oxide concentration. 

The rate law governs the conversion the nitric oxide and 
oxygen nitrogen 


= k [NO]? [O.] 


The time required for conversion half the nitric oxide may 
calculated from the equation 


since the oxygen concentration the tunnel was found 
essentially constant. 

The concentrations total oxides nitrogen relation 
nitrogen dioxide found the tunnel atmosphere reflect the 
relatively slow rate oxidation nitric oxide 
dioxide. This result was anticipated, since the primary oxidation 
product the atmospheric nitrogen supplied for combustion 
with the diesel fuel oil nitric oxide. Hence, was expected 
that nitric oxide concentration the tunnel would tend 
greatly exceed those nitrogen dioxide. Table shows typical 
concentration levels, from which clear that several hours 
would required for conversion appreciable quantities 
nitric oxide the much more toxic nitrogen dioxide. 

Under normal operation, most the gases are swept out 
the tunnel the ventilation caused the piston effect 
passing trains and view the low concentrations nitrogen 
dioxide, was considered that concentrations oxides 
nitrogen presented serious hazard. 


(d) Smoke and Particulate Matter 


Smoke and particulate matter occurred relatively high 
concentrations the tunnel environment, especially during 
periods heavy train activity. Smoke concentrations deter- 
mined hourly basis, attained values considerably higher 
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Date Samples 


Time 


4p.m. 
4p.m. 


Highest 


values exceeding 
Concentration 


Mean M.A.C. p.p.m. 
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Exhaust Constituent Mean 


Total Oxides 


4p.m. 
p.m. p.m. 3.6 


Highest 


Values 
Concentration 


Exceeding M.A.C. 


M.A.C. 
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than those usually encountered both urban and industrial 
The mean hourly smoke concentrations for 
hr. periods varied from 5.3 10.1 COH units per 1,000 
linear air, with maximum hourly concentration 22.4 
linear ft. One COH unit defined that quantity 
light scattering solids producing optical density 0.01 
when measured light transmission. Parallel observations 
sampling station outside the tunnel, Station yielded values 
0.4 3.0 COH/1,000 linear ft. whereas Station inside 
the tunnel the values ranged from 2.2 COH/1,000 linear 
fr. (Figure 

variation concentration (over hr. test periods) from 
0.9 2.3 mg./cu. was found for the particulate matter 
collected the high volume filtration samplers. Typical values 
for particulate matter collected during the test days are 
shown Table 

Smoke and particulate matter represent non-specific air 
pollutants that are objectional above certain minimum levels 
depending upon effects visibility and possible respiratory 
effects. relatively high values obtained for smoke and 
particulate factors considered the subsequent 
decision install forced ventilation the tunnel. 


(e) Air Flow the Tunnel 


The displacement part the tunnel atmosphere moving 
trains results piston action which pushes considerable 
volume air ahead the train while corresponding suction 
effect draws fresh air into the tunnel the inbound portal. 
Smoke and exhaust gases the form trailing wake 
“slug” are drawn along with the train. (Figure 

Measured average air displacement per train was 1.54 
cu. ft., which somewhat less than the calculated volume 
the tunnel 1.87 ft. Therefore, since, the average, 
less than one complete air change occurred the tunnel 
result the passage train, low residual concentrations 
exhaust components were found the tunnel atmosphere. These 
were determined from the samples, that samples 
collected prior entrance trains into the tunnel. 
residual levels toxic gases the tunnel may also computed 


when the residual gas concentration, the original gas 
concentration, the coefficient ventilation and the 
persistence time. The number air displacements necessary 

For example, where initial concentration total oxides 
nitrogen was 47.2 p.p.m. the tunnel’s midpoint, and 


| COH REPRESENTS CLEAR ATMOSPHERE 
10 COHS IS DISAGREEABLY SMOKY 


wo 


Figure 4—Hourly smoke concentrations Stations and 


October 23, 1957. 
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TABLE 
CONCENTRATIONS PARTICULATE MATTER MILLIGRAMS 


METRE CALCULATED FROM 
VOLUME FILTER WEIGHINGS 
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maximum air displacement was 2.6 10° cu. ft., calculation 
residual concentration after fresh air was drawn 
train’s passage showed the residual concentration 11.9 
Where air displacement was minimal, namely 1.0 
cu. ft., total oxides nitrogen were recorded concentration 
25.2 p.p.m. Calculation reduces this residual level 
concentration 14.8 p.p.m. 


Ventilation the Tunnel 


Four axial flow fans were installed, two each portal 
shown Figure with one fan each side the portal. 


Each these fans has output 57,000 c.f.m. against free 
air 47,000 c.f.m. against resistance in. water. When 
the fans are reversed they can withdraw 60% these quantities. 
Each fan powered h.p., 1,800 r.p.m., cycle, phase 
550 motor. The motor the totally enclosed type directly 
coupled the fan within the fan housing. 


The fans are in. diameter the outside the housing. 
Each fan outlet extended into the tunnel approximately 
means air duct, approximately sq. ft. cross-sec- 
tional area. The fans are located just outside the portals and the 
ducts are fastened the tunnel walls maintain full estab- 
lished clearances for the passage locomotives and trains. 


The operation the fans controlled photocells located 
few feet outside each portal. The passage train starts 
all fans blowing air through the tunnel the direction train 
movement. time delay relay permits the fans continue 
operating for pre-set interval after the train has cleared the 
outbound portal. This interval adjusted from minutes 
and limited only the practical interval dispatching the 
next train through the tunnel. Operation the ventilation 
system supervised remote indication and control the 


Figure 5—St. Clair tunnel—smoke the form “Slug” 
drawn along with the train. 
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Figure 6—Ventilating fan. There are two axial flow fans 
installed each portal. 


dispatchers’ tower the Sarnia tunnel summit, and malfunction 
the system can therefore detected immediately. 

The ventilating system induces air velocity approx- 
imately 3.5 m.p.h. the tunnel and designed supplement 
the ventilation produced the piston action trains. This 
system somewhat novel compared the high velocity 
ventilation which has been installed several railroad tunnels. 
Employment this type system was possible due the 
peculiar shape and construction the tunnel, described 
which results considerable amount ventilation 
induced train movement and the absence natural air flow 
the tunnel result variable winds and other weather 
conditions. 

The combined capacity the four fans that provide ventilat- 
ing air the tunnel equal volume air 
permitting complete change air the tunnel every 
minutes, this ventilating effect being addition that induced 
train passage. 

Extensive tests conducted the tunnel subsequent the 
installation the ventilating system have shown that concentra- 
tions smoke and oxides nitrogen are reduced approx- 
imately 30% with corresponding reduction concentrations 
the other constituents diesel exhaust gas. 

These reductions are sufficient bring the concentrations 
all components the exhaust fumes well below acceptable 
levels and provide adequate safety margin for all conditions 
operations. 

Following the installation the forced ventilating system 
and tests evaluate its performance, full dieselized operation 
the tunnel was authorized the fall 1958. 
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